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; Abstract 

•rH , 

X ■ 

; i ■ Valuable data on quarkonia (the bound states of a heavy quark Q = c,b and 

the corresponding antiquark) have recently been provided by a variety of sources, 
mainly e + e~ collisions, but also hadronic interactions. This permits a thorough 
updating of the experimental and theoretical status of electromagnetic and strong 
transitions in quarkonia. We discuss QQ transitions to other QQ states, with some 
reference to processes involving QQ annihilation. 
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I. INTRODUCTION 



Quarkonium spectroscopy has celebrated a great resurgence in the past few years thanks 
to a wealth of new information, primarily from electron-positron colliders, but also from 
hadronic interactions. Transitions between quarkonium states shed light on aspects of quan- 
tum chromodynamics (QCD), the theory of the strong interactions, in both the perturbative 
and the non-perturbative regimes. In the present article we review the new information on 
these states and their tr ansitions and indicate theore t ical implications , upda ti ng earlier dis- 



cussions such as those in 



(2004): 



Barnes and Godfrey! (120 04); 



Godfrey and Rosnerl ( 



2001a 



Brambilla 



Kwong et all ( 119881 . Il987f ) 



et al 



(120041): 



Eichten et al. 



Kwong and Rosner 



( 119881 ) (which may be consulted for explicit formulae). 

We shall deal with states composed of a heavy quark Q = c or b and the corresponding 
antiquark Q. We shall discuss QQ transitions primarily to other QQ states, with some ref- 



erence to processes invol ving QQ annihilation, and w i 



treat e d, for example , in 



( 120041 ): 



Eichten et al 



Barnes and Godfrey! (120041 ) : 



(12004 . 



1 largely bypass dec a ys to open flavor 



Barnes et al. 



(|20Q5| ); 



Brambilla et al. 



20061 )). 



A brief overview of the data on the charmonium and T systems is provided in Section [TT1 
We then review theoretical underpinnings in Section II III discussing quarks and potential 
models, lattice gauge theory approaches, perturbative QCD and decays involving gluons, 
and hadronic transitions of the form QQ — > (QQ)'+ (light hadrons). Section [TV] is devoted 
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\j/ (4S) or hybrid 

*» DD 

V(3S) 




- j pc = - + i~ i + - o ++ i ++ 2 

L = 1 1 1 1 



FIG. 1 Known charmonium states and candidates, with selected decay modes and transitions. Red 
(dark) arrows denote recent observations. 

to charmonium and Section |V] to the bb levels and includes a brief mention of interpolation 
to the be system. Section [VI] summarizes. 

II. OVERVIEW OF QUARKONIUM LEVELS 

Since the discovery of the J/ip more than thirty years ago, information on quarkonium 
levels has grown to the point that more is known about the cc and bb systems than about 
their namesake positronium, the bound state of an electron and a positron. The present 
status of charmonium (cc) levels is shown in Fig. [U while that of bottomonium (bb) levels 
is shown in Fig. [2J The best-established states are summarized in Tables [U and HT1 

The levels are labeled by S, P, D, corresponding to relative orbital angular momentum 
L = 0, 1, 2 between quark and antiquark. (No candidates for L > 3 states have been 
seen yet.) The spin of the quark and antiquark can couple to either S = (spin- singlet) 
or S = 1 (spin-triplet) states. The parity of a quark-antiquark state with orbital angular 
momentum L is P = (— 1) L+1 ; the charge-conjugation eigenvalue is C = (— 1) L+S . Values 



5 




FIG. 2 Transitions among bb levels. There are also numerous electric dipole transitions S «-> P «-> D 
(not shown). Red (dark) arrows denote objects of recent searches. 

of J PC are shown at the bottom of each figure. States are often denoted by 2S+1 [L]j, with 
[L] = S, P, D, . . .. Thus, L = states can be 1 S'o or 3 Si; L = 1 states can be 1 P\ or 3 -Po,i,2! 
L = 2 states can be 1 D 2 or 3 -Di 5 2,3, and so on. The radial quantum number is denoted by n. 



III. THEORETICAL UNDERPINNINGS 



A. Quarks and potential models 



An approximate picture of quarkonium states may be obtained by describing them as 
bound by an interquark force whose short-distance behavior is approximately Coulombic 
(with an appropriate logarithmic modification of coupling strength to account for asymptotic 
freedom) and whose long-distance be havior is linear to account for quark co nfinement. An 



example of this approach is found in 



may be found in 



Appelauist et al. 



Eichten et al. 



1978) 



Grossc and Martin 



Quigg and Rosnerl (119791 ). iRadford and Repkol (120071 ) presents more recent results. 



197ol 197 j 197 j . 198(4): early 



(1198 



3); 



reviews 



Novikov et al. 



tog); 
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TABLE I Observed charmonium states. All numbers are quoted from lYao et al\ (|2006l ). More 
recent information is included in the text, where available. 



Quantum numbers 


Name 


Mass 


Width 


n 


L 


J PC 


n 25+1 Lj 




MeV 


MeV a 


1 





o-+ 






2980.4 ± 1.2 


25.5 ±3.4 


1 





1— 






3096.916 ±0.011 


93.4 ± 2.1 keV 


1 


1 


0++ 


l 3 Po 


Xoo(lP) 


3414.76 ± 0.35 


10.4 ±0.7 


1 


1 


1++ 


l 3 Pi 


Xcl(lP) 


3510.66 ± 0.07 


0.89 ±0.05 


1 


1 


2++ 


1 3 P 2 


Xc2(lP) 


3556.20 ± 0.09 


2.06 ±0.12 


1 


1 


1+- 


I'Pi 




3525.93 ± 0.27 


< 1 


1 


2 


1— 


l 3 D 1 


^(3770) 


3771.1 ±2.4 


23.0 ±2.7 


2 





o-+ 


2 1 S 


»? c (25) 


3638 ± 4 


14 ±7 


2 





1— 




^(25) 


3686.093 ± 0.034 


337 ± 13 keV 


2 


1 


2++ 


2 3 P 2 


Xc2(2P) 


3929 ± 5 


29 ± 10 



a Unless noted otherwise. 



1. Validity of nonrelativistic description 

In order to estimate whether a nonrelativistic (NR) quarkonium description makes sense, 
"cartoon" versions of cc and bb spectra may be constructed by noting that the level spacings 
are remarkably similar in the two cases. They wo uld be exactly equal if t he interquark 
potential were of the form V(r) = C\og(r/r ) (see Quigg and Rosner ( 1977 7 ! ) ) . which may 
be regarded as a phenomenological interpolation between the short-distance ~ — 1/r and 
long-distance ~ r behaviors expected from QCD. In such a potential the expectation value 
of the kinetic energy (T) = (r/2){dV/dr) is just C/2 ~ 0.37 GeV with C = 0.733 as found 



in 



Quigg and Rosnerl (11979 ). Since (T) = 2 • (l/2)mQ(t> 2 ), one has (v 2 ) ~ 0.24 for a charmed 



quark of mass m c ~ 1.5 GeV/c 2 (roughly half the J/ip mass), and (v 2 ) ~ 0.08 for a b quark 
of mass nib ~ 4.9 GeV/c 2 (roughly half the T(1S) mass). Thus a nonrelativistic description 
for charmonium is quite crude, whereas it is substantially better for bb states. 
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I 1 J — I 

TABLE II Observed bottomonium states. All numbers are quoted from lYao et al\ (|2006l ). More 
recent information is included in the text, where available. 



Quantum numbers 


Name 


Mass 


Width 


n 


L 


J PC 


n 2S+1 Lj 




MeV 




1 





1~ 


1 3 S 1 


T(LS) 


9460.30 ± 0.26 


54.02 ± 1.25 keV 


1 


1 


0++ 


l 3 Po 


Xbo(lP) 


9859.44 ± 0.52 


unknown 


1 


1 


1++ 


l 3 Pi 


XbiO-P) 


9892.78 ± 0.40 


unknown 


1 


1 


2++ 


1 3 P 2 


XiaO-P) 


9912.21 ±0.40 


unknown 


1 


2 


2™ 


l 3 P>j a 


T(ID) 


10161.1 ± 1.7 


unknown 


2 





\— 


2 3 5i 


T(25) 


10023.26 ± 0.31 


31.98 ± 2.63 keV 


2 


1 


0++ 


2 3 P 


X6o(2P) 


10232.5 ± 0.6 


unknown 


2 


1 


1++ 


2 3 Pi 


X6i(2P) 


10255.46 ± 0.55 


unknown 


2 


1 


2++ 


2 3 P 2 


X6 2 (2P) 


10268.65 ± 0.55 


unknown 


3 





1~ 


3 3 S 1 


T(35) 


10355.2 ± 0.5 


20.32 ± 1.85 keV 


4 





1~ 


4 3 5i 


T(45) 


10579.4 ± 1.2 


20.5 ± 2.5 MeV 



"Probably all or mostly J = 2. 
2. Role of leptonic partial widths: |\&(0)| 2 

The partial widths for 3 Si states to decay to a lepton pair through a virtual photon are 
a probe of the squares l^lO)) 2 of the rela tive n 3 Si wave functions at the origin through the 
relation (IVan Royen and WeisskopJ . 119671 ) 



r(n 3 Si -> e + e 



167ra 2 eg|^ r 



x 1 



16«5 

3vr 



(1) 



where cq = 2/3 or —1/3 is the quark charge , M n i s the mass of the n 3 Si state, and the 
last term is a QCD correction (IKwong et all Il988l ). Thus leptonic partial widths are a 
probe of the compactness of the quarkonium system, and provide important information 
complementary to level spacings. Indeed, for the phenomenologically adequate potential 
V(r) = Clog(r/r ), a change in the quark mass tuq can be compensated by a change in 
tq without affecting quarkonium mass predictions (tq can be viewed as setting the overall 
energy scale), whereas a larger quark mass will lead to a spatially more compact bound state 
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and hence to an increased value of |\I/(0)| 2 for each state. A more general form is the power- 
law potential, V(r) ~ sgn(u)r u , which approaches the logarithmic potential in the limit of 
v — > 0. One can show that in the power-law potential lengths sc ale as m^ 1 ^ 2 ^ an d henc e 



\]/(0)| 2 scales as m 3 Q^ 2+u \ or ~ m^m^ ,uiq for v = —1,0,1 (IQuigg and Rosnerl . Il979l ) 



3/2 



(In charmonium and bottomonium the gro und states have sizes of about 0.4-0.5 fm and 0.2 
fm, respectively ( Quigg and Rosnerl . 1981 ).) Thus the effective quark mass in a potential 
description is constrained by measured leptonic widths. One can expect that in descriptions 
such as lattice gauge theories, to be discussed in Section IIII.B| similar constraints will hold. 

The scaling of leptonic widths from the charmonium to the bottomonium family can 
be roughly estimated using the above discussion, assuming an effective power v ~ 0. In 
that case the leptonic width for each n scales as T ee (nS) oc eql^ (0)\ 2 / uiq oc e^/rriQ 2 . As 
the QCD correction in Eq. (1) is appreciable [as are relativistic corrections, particularly for 
charmonium], this is only an approximate rule. 

The important role of leptonic widths is particularly ev ident in constructi o ns of 



the interquark pq t entia 



Quigg and Rosner . 



1981 



based on inverse-scat t ering methods 



Schonfeld et al. 



198 



0; 



Thacker et al. 



( Kwong and Rosnerl . 



1978a 



1986; 



bj). The reduced ra- 



dial wave functions u n s(r) = r^/ n s(r) on the interval < r < oo for an S'-wave Schrodinger 
equation with central potential V(r) may be regarded as the odd-parity levels (since they 
must vanish at r = 0) in a symmetric potential V(—r) = V(r) on the interval — oo < r < oo. 
The even-parity levels [with u(0) ^ 0] do not correspond to bound states but, rather, equiva- 
lent information is provided by the leptonic widths of the nS levels, which gives the quantities 



I* 



= |w^g(0)|. Thus, if QCD and relativistic corrections can be brought under control, 
leptonic widths of the S'-wave levels are every bit as crucial as their masses. 

A recent prediction of the l eptonic width ratio r ee [Y(2S)]/r ee [T(lS)] = 0.43 ± 0.05 in 
lattice QCD ( Gray et"q/] . 2005 ) raises the question of what constitutes useful measurement 
and prediction precisions, both for ratios and for absolute leptonic wid ths. (For compariso n, 
the CLEO Collaboration has measured this ratio to be 0.457 ± 0.006 (jRosner et all 120061 ).) 
Potential models have little trouble in predicting ratios T ee (n'S)/T ee (nS) to an accuracy of 
a few percent, and one would thus hope for lattice approaches eventually to be capable of 
similar accuracy. Much more uncertainty is encountered by potential modes in predicting 
absolute leptonic widths as a re sult of QCD and re l ativist ic corrections (see, for example, the 
inverse-scattering approach of iQuigg and Rosnerl (1l98ll )). Measurements with better than 
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a few percent accuracy, such as those in iRosner et al\ (120061 ) and others to be discussed 
presently, thus outstrip present theoretical capabilities. 



3. Spin- dependent interactions 



Hyperfine and fine-structure splitting s in quarkonium are s e nsitive to the 



structure of the interquark interaction (lAppelquist et al 



Kwong et al 



1978 



Brambilla et al. 



jorentz 



1987 



Novikov et al 



2004 : 



19781 ). One may regard the effective potential V(r) as 



the sum of Lorentz vector Vy and Lorentz scalar V$ contributions. The spin-spin interaction 
is due entirely to the Lorentz vector: 



V ss (r) 



(2) 



where <jq and <jq are Pauli matrices acting on the spins of the quark and antiquark, re- 
spectively. For a Coulomb-like potential ~ — 1/r the Laplacian is proportional to S 3 (r), so 
that Vss{r) contributes to hyperfine splittings only for S waves, whose wave functions are 
non-zero at the origin. In QCD the coupling constant undergoes slow (logarithmic) varia- 
tion with distance, leading to small non-zero contributions to hyperfine splittings for L > 
states. Relativistic corrections also result in small non-zero contributions to these splittings. 
Both spin-orbit and tensor forces affect states with L > 0. The spin-orbit potential is 



V LS (r) 



L-S fdV v dV f 



2m 2 r\ r \ dr 



dr 



(3) 



where L is the relative orbita 



spin. The tensor potential is (IMessiah 



angular momentum of Q and Q, while S is the total quark 



1999 



Radford and Repko 



20071 ) 



Vr(r) 



s 7 



12mj 



1 dV v d 2 V v 



r dr 



dr 2 



(4) 



with St = 2[3(S -f)(S -f) — S 2 } (where S = Sq + Sq is the total spin operator and f is a unit 
vector) has non-zero expectation values only for L > [e.g., —4,2, —2/5 for 3 -Po,i,2 states]. 

B. QCD on the lattice 

At momentum scales less than about 2 GeV/c (distance scales greater than about 0.1 fm) 
the QCD coupling constant o>s{Q 2 ) becomes large enough that perturbation theory cannot 
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be used. The value a s {m 2 T ) = 0.345 ± 0.010 (iBethkd . 120071 : bavier et all 120071 : iKluthl . |200J) 
is just about at the limit of usefulness of perturbation theory, and as(Q 2 ) increases rapidly 
below this scale. One must resort to non-perturbative methods to describe long-distance 
hadronic interactions. 

If space-time is discretized, one can overcome the dependence in QCD on perturbation 
theory. Quark confinement is established using this lattice gauge theory approach. For low- 
lying heavy quarkonium states, below the threshold for Zweig-allowed decay to open heavy 
flavor mesons, an accurate description of the spectrum can be obtained, once one takes 
account of the degrees of fr eedom associated with the production of pairs of light (u, d, s) 
quarks (IDavies et all 120041 ). For exam ple, recent lattice ca lculations of the spin-splitting 



20071 ). while the experimental value 



between J/ ip and rj c yield 1 1 1 ± 5 MeV (IFollana et al 
is 117.1 ± 1.2 MeV. 

Above threshold, the situation is more challenging: (1) Heavy quarkonium states have 
more typical-size hadronic widths; (2) Such states are usually not the ground state for a given 
set of quantum numbers; and (3) These resonances are embedded in a multibody continuum. 
In the lattice approach, information is extracted from Euclidean correlation functions. This 
makes dealing with excited-state resonances in a multibody continuum particularly difficult 



teulava et al. 



20071 ). 



Lattice QCD also provides a theoretical underpinning for the phenomenological potential 
model approach. The well-measured static ener gy be t ween a heavy quark-antiquark pair 
justifies the form of the nonrelativistic potential (jBalil . 1200 ll ). Recently, high- accura c y lat 



tice calculations o 



Koma and Komal . 



the spin-dependent potentials have also been made (iKoma et al. 



2006 



20071 ). This approach allows the direct determination of the spin-orbit, 



spin-spin and tensor potentials as well. At present, these spin-dependent potential calcula- 
tions have not yet included the effects of light quark loops. 



C. Electromagnetic transitions 

The theory of electromagnetic (EM) transitions between quarkonium states is straightfor- 
ward, with terminology and techniques familiar from the study of EM transitions in atomic 
and nuclear systems. Although electromagnetic transition amplitudes can be computed from 
first principles in lattice QCD, these calculations are in their infancy. At the present time, 
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only potential model approaches provide the detailed predictions that can be compared to 
experimental results. In this approach, the spatial dependence of EM transition amplitudes 
reduces to functions of quark position and momentum between the initial and final state 
wave functions. Expanding the matrix elements in powers of photon momentum generates 
the electric and magnetic multipole moments and is also an expansion in powers of veloc- 
ity. The leading order transition amplitudes are electric dipole (El) and m agnetic dipole 



(Ml) . In what follows we shall take m c = 1.5 GeV/c 2 and m;, = 4.9 GeV/c 2 (IKwong et al. 



19881 ). which are considered "constituent-quark" values, appropriate to the non-perturbative 



regime found in charmonium and bottomonium. 



1. Magnetic dipole transitions 



Magnetic dipole transitions flip the quark spin, so their amplitudes are proportional to 
the quark magnetic moment and therefore inversely proportional to the constituent quark 
mass. At leading order the magnetic dipole (Ml) amplitudes between 5-wave states are 
independent of the potential model: The orthogonality of states guarantees, in the limit of 
zero recoil, that the spatial overlap is one for states within the same multiplet and zero for 
transitions between multiplets which have different radial quantum numbers. 

Including relativistic corrections due to spin dependence in the Hamiltonian spoils this 
simple scenario and induces a small overlap between states with different radial quantum 
numbers. Such n ^ n' transi tions are referred to as "hindered" . In cluding finite size correc- 



tions the rates are given by (jEichten et al. 



1975 



1976 . 



1978 



1980) 



Aae^(2J f + l)\(f\Ukr/2)^/3mi 



(5) 



T{n 3 S 1 -W 1 S + 7 ) 
r(n l S -+ n^Sx + 7) 

where eg = 2/3 or —1/3 is the quark charge, k is the photon energy, jo(x) = sinx/x, and 
rriQ is the quark mass. The only Ml transitions between quarkonia states so far observed 
occur in charmonium, but the corresponding transitions in bb systems are the objects of 
current searches. For small k, jo(kr/2) — > 1, so that transitions with n' = n have favored 
matrix elements, though the corresponding partial decay widths are suppressed by smaller 
k 3 factors. 

Numerous papers have studied these Ml transitions including full relativistic corrections 
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(Ebei 


t et al. 


2003a 




Godfrey and Isgur 


, 1 


985 


; 


Glodfrev and 


losner 


1984: 


Lahde. 


20031: 


Zambetakis and Bvers. 


1983 


Zhang et al.. 


1991) 



2001b 



Grotch et al. 



19911 ). They depend explic- 



itly on the Lorentz structure of the nonrelativistic potential. Several sources of uncertainty 
make Ml transitions particularly difficult to calculate. In addition to issues of relativistic 
corrections and what are known as "exchange currents," there is the possibility of an anoma- 
lous magnetic moment of the quark (kq). Furthermore, the leading-order results depend 
explicitly on the constituent quark masses, and corrections depend on the Lorentz structure 
of the potential. 



2. Electric dipole transitions 



The partial widths for electric dipole (El) transi t ions between states 2S+1 3Lij. and 



2S+1 Lf j^ are given by (jEichten et al. 



1975 



1976 



19781 . Il980h 



T(n 2S+1 L 



- h 



n 



/25+1 



<fj, + 7) 



4aegk 3 



(2J f + l)S if \{f\r\i)f 



(6) 



The statistical factor <Sj/ is 

S if = S fi = max(Lj, Lj 



Ji 1 Jf 
Lf S Li 

For transitions between spin-triplet S"-wave and P-wave states, Sif — , r 
P <-> D transitions, which have also been observed both in charmonium and in the hh 



(7) 

q. Expressions for 



system, are given, for example, in lKwong and Rosnerl (119881 ) . 



The leading corrections f o r elect r ic dipole corrections have 



ber of authors (Ebert et al 



Grotch et al 



1984; 



Mc Clary and Byers 



2003a 



Feinberg and Sucherl. 



Grotch and Sebastian 



198 



3: 



1982 



Moxhay and Rosner 



j een considered by a num 



1975 



Godfrey and Isgur 



Kang and Sucherl. 



1983 



Suchei 



19781 : 



1978 



Lahd 



1985 



2003; 



Zambetakis and Bvera . 



Grotch et al. 



198J). 



19831 ). A general form was derived by Grotch, Owen and Sebastian ([( 
There are three main types of corrections: relativistic modification of the nonrelativistic wave 
functions, relativistic modification of the electromagnetic transition operator, and finite-size 
corrections. In addition to these there are additional corrections arising from the quark 
anomalous magnetic moment. For the 3 Pj <-> 3 Si transitions in which we are primarily 
interested, the dominant relativistic corrections arise from modifications of the wavefunc- 



tions and are included by the quarkonium analog of Siegert's theorem (IMcClary and Byersl . 
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1983 



Siegertl . 119371 ) . We will find that differences in theoretical assumptions of the various 



potential models make it difficult to draw sharp conclusions from the level of agreement 



of a particular model with experimental data. However, there is usua 



varia tion in the NR predictions if the models are fit to the same states (IKwong and Rosnerl . 



ly very little model 



19881 ). The only exceptions are transitions where the dipole matrix element exhibits large 
dynamical cancellations, for instance when higher radial excitations are involved which have 
nodes in their wavef unctions. 



3. Higher multipole contributions in charmonium 



Magnetic quadrupole (M2) amplitudes are higher order in v 2 /c 2 . They are of inter- 
est because they provid e an indirect measure of the charmed quark's magnetic moment 



( jKarl et al. 



1976 



1980j) and are sensitive to D-w&ve admixtures in S"-wave states, pro- 



viding another means of studying the l 3 £>i - 2 Z S X mixing in the ^(3770) — ip(2S) states 



(jGodfrey et al 



1986 



Sebastian et al 



19921 ). They affect angular distributions in decays 
such as 4>(2S) — > Xcj + 1 and Xcj J/ip + y and become experimentally accessible through 
interference with the dominant El amplitudes. 

The Xcj iJ/i* or ^'(25') — > ^Xcj decays may be described by the respective helicity 
amplitudes A\ or A' x , in which A labels the projection of the spin of the Xcj parallel (for A\) 
or antiparallel (for A' x ) to the photon, which is assumed to have helicity +1. The radiative 
widths are given in terms of these amplitudes by 



IW2S) - 

r(xcj 



IXcj) 



j/4>) 



£14 



A>0 

El 



2J 



(8) 
(9) 



A>0 



In terms of a parameter e = ^£? 7 /(4m c ), where ^ = — 1 for ip(2S) — > 7Xcj and ^ = +1 
for Xcj — > yJ/i>i the predicted helicity amplitudes A\ or A' x are in the relative proportions 
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(Karl et al. 


> 


1976 




1980) 



X c2 : A 2 = VQ{1 + e(l + K c )] (10) 

A 1 = v/3[l - e(l + k c )\ (11) 

A = [1 - 3e(l + k c )] (12) 

Xci : A = >/3[l + e(l + K c )] (13) 

A = V3[l - e(l + k c )] (14) 

Xc0 : A = V2[l - 2e(l + K c )] . (15) 



Here an overall El amplitude has been factored out, and k c is the charmed quark's anomalous 
magnetic moment. 



D. Perturbative QCD and decays involving gluons 



Many quarkonium decays proceed through annihilation of QQ into final states consisting 
of gluons and possibly photons and light-quark pair s. Expressions for partial w idths of color 



Petrelli et al. 



single t QQ systems are given in lKwong et al\ (119881 ). and have been updated in 
( 119981 ) . In that work, annihilation rates are also given for the color-octet component of 
the QQ system, which appears necessary for successful description of QQ production in 
hadronic interactions. We shall confine our discussion to the effects of the color-singlet QQ 
component in decays. Discrepancies between theory and experiment can be ascribed in part 
to neglected relativistic effects (particularly in charmonium) and in part to the neglected 
color-octet component. 



E. Hadronic transitions [QQ — * (QQ)'+ (light hadrons)] 

A number of transitions from one QQ state to another occur with the emission of light 
hadrons. So far, the observed transitions in charmonium include ip{2S) — > J/ip7T + 7T~, 
ip(2S) -> J/^ttV, 4){2S) J/ipT], ip(2S) -> J/ipir , and ip(2S) -> h c ir°. In addi- 
tion, above charm threshold a state X(3872) decays to J/ipir + 7t~, and a state F(3940) 
decays to J/ipu. The observed transitions in the bb system include T(2S) — > T(lS)nir, 
T(3S) -> T(15,25)7T7r, x(2P)bi, 2 -> T(lS)u, and x(2P) bJ -> Xwtttt. Many of these tran- 
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sitions have been observed only in the past few years (see later sections for experimental 
data). 



The theoretical des cription of had r onic 



transitions uses a mu 



emission deye l oped in 



Peskinl (119791 ): 



Voloshin 



ranot et all 



lazsl); 



ti pole e xpansion for gluq n 



1979|); Bhanot and Peskinl (Il979h : iGottfriedl fll978h : 



Yanl (119801 ). Formally, it resembles the usual multipole ex- 



pansion for photonic transitions discussed in Section IIII.CI The interaction for color electric 
and magnetic emission from a heavy quark is given by 



H I = y"rf 3 xg t (x)t a [x-E a (x) + ( r-B a (x)]g(x) + ... , 



(16) 



where t £ 



8) is a generator of color SU(3) , and the ( Q)Q and E, B are dressed 



(anti) quarks and color electric and mag netic fields (jYarj , ll98pL As usual, the multipole 
expansion arises from expanding the color-electric and color-magnetic fields about their 
values at the center of mass of the initial quarkonium state. However, unlike EM transitions, 
a single interaction of Hj changes a color singlet QQ initial state (i) into some color octet 
QQ state. Therefore, a second interaction Hj is required to return to a color singlet QQ final 
state (/). In the overall pro cess at least two gluons are emitted. Assuming factorization 



for the quarkonium systems (IKuang and Yanl . Il98ll ). the full transition amplitude can be 
expressed as a product of two subamplitudes: One that acts on the quarkonium system to 
produce the multipole transition and a second that creates the final light hadrons (H) from 
the action of the gluonic operators on the vacuum state - 



In non-relativist ic QCD (NRQCD) (IBodwin et al 



Luke and Manohar 



1995 



Caswell and Lepagd . [1986; 



19971 ). the strength of the various interactions can be ordered in powers 



of the heavy quark velocity v. The leading behavior comes from t wo color-electric (El) glu on 



emissions. This amplitude can be written in the factorized form (IKuang and Yan 



19811 ): 



E< ^t1O)<O|r^|/) (0|EiE£|g) 



o 



Ej — E, 



(17) 



The sum runs over allowed QQ octet intermedi ate states O. Phenqmenol ogical models 
(e.g. the Buchmuller-Tye vibrating string model ( Buchmuller and Tye . 198ol )) are used to 
estimate this quarkonium overlap amplitude. The quantum numbers of the initial and final 
quarkonium states determine which terms in the multipole expansion may contribute. For 
the light hadron amplitude the the states allowed are determined by the overall symmetries. 
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In transitions between various 3 Si quarkonium states the leading term in the multipole ex- 
pansion has two color-electric (El) interactions. The lowest-mass light hadron state allowed 
is a two-pion state with either an S- or D-w&ve relative angular momentu m. The form of the 
light hadron amplitude is determined by chiral symmetry considerations (IBrown and Cahnl . 
19751 ): 



(0|EiE^|7r(A; 1 )7r(A; 2 )) = 5^0^% ■ k 2 + c 2 (k{k k 2 + k\k\ - ■ k 2 )}. (18) 

The two unknowns (ci, c 2 ) are the coefficients of the S-wave and D-wave two-pion systems. 
Thei r values are det ermined from experiment. Additional terms can arise in higher orders 
in v jvoloshinl . hoOffl . 

Hadronic transitions which can flip the heavy quark spins first occur in amplitudes with 
one color-electric (El) and one color-magnetic (Ml) interaction. These transitions are sup- 
pressed by an additional power of v relative to the purely electric transitions. Transitions 



involving two color-magnetic interactions (Ml) are suppressed by an a dditiona 



Many detailed pred i ctions for hadron i c transition r ates can be found in 



Kuang et all f|l988l ): iKuang and Yanl (Il98ll . Il990l ): IVoloshinl (11986 



2003 



Kuans! (12002 



power of y . 



20061 ). 



IV. CHARMONIUM 



In what follows we shall quote masses and partial widths from lYao et al\ (120061 ) unless 
otherwise noted. The masses are used to calculate photon transition energies. We shall 
use an electromagnetic coupling constant a = 1/137 in all cases. For gluon emission in QQ 
annihilation we shall use a momentum-dependent strong coupling constant cxs{Q 2 ) evaluated 
at Q 2 = rrin. The QCD cor rections to the decay widths we quote are performe d for this scale 



choic e (IKwong et al. 



19881 ). Typical values are asijn 2 c ) ~ 0.3, otsijn 2 ) ~ 0.2 (|Kwong et al 



1988). A different scale choice would lead to different 0(as) corrections ( jBrodsky et al 
1983k 



A. The J /if) 



The J/jj w a s the first charmonium state discovered, in 1974 (jAubert et al. 



1974 : 



Augustin et all 1 1 9 T4f ) . It is the lowest 3 Si cc state and thus can couple directly to vir- 
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tual photons produced in e + e colli sions. The most precis e mass determination to date 



comes from the KEDR Collaboration (lAulchenko et al 



2003|), m(J/ip) = 3096.917±0.010± 



0.007 MeV, a relative uncertainty of 4 x 10 6 . 

The J/ip intrinsic width originally w as determined in directly. The history of these mea- 



surements shows values below 70 MeV (IBai et al. 
suring the excitation curve in pp — > e + e~ ([Armstrong et al. 



1995 ). A di r ect d etermination by mea- 



)3) 



1993) was the first to re- 



sult in a substantially higher value, albeit still with considerable statistical uncertainty: 
F(J/if)) = 99 ± 12 ± 6keV. Recent indi r ect measurements, r e sulting in uncertainties of 



3 — 4 keV, were carried out (lAdams et all l2006bl ; lAubert et all l2004al ) using the radiative 
return process e + e~ — > r je + e~ — > 7J/V — > 7(/U + /i~). The experimental observable is the 
radiative cross-section, a convolution of the photon emission probability and the J/ ip Breit- 
Wigner resonance shape. It is calculable and proportional to the coupling of the J/ip to 
the annihilating e + e~ pair and the J/ip decay branching fraction, T ee x B(J/if> — > fi + fi~). 
Interference with the QED process e + e~ — > 7/i + /i~ introduces an asymmetry around the 
J/i/j peak in m(/z + /i~) and must be taken into account. B(J/ip — > is known well, 

hence the product gives access to Y ee and, toget her with B j Jjjb -^ /i + /i~), to r to t- The 



current world average is T(J 'fi]j) = 93.4 ± 2.1 keV fjYao et all 120061 ) 



The largest data sample now consists of 58 million J/ijj collected by the BES-II Collabora- 
tion. Decays from the ip(2S) state, in particular ip(2S) — > 7T + 7T~ J/tp — > 7r + 7r _ hadrons, offer 
a very clean avenue to study J/ip final states, yielding one iT + 7t~J/tp event per three ip(2S) 
produced. Experimentally, this can be handled by requiring a tt + tt~ pair recoiling against 
a system of m(J/ip), without further identification of the J/ip decay products. This path 
also eliminates contamination of the sample by continuum production of a final state under 
study, e + e~ — > 7* — > hadrons. Other J ftp production mechanisms include pp collisions and 
radiative return from e + e~ collisions with center-of-mass energy > m(J/%p). Many decays 
of J/i/j to specific states of light hadrons provide valuable information on light-hadron spec- 
troscopy. Here we shall be concerned primarily with its decay to the ?7 c (l 1 5'o), the lightest 
charmonium state of all; its annihilation into lepton pairs; and its annihilation into three 
gluons, two gluons and a photon, and three photons. 
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1. J/ip -> 7% 



The rate predicted for the process J/tj) — > ^r\ c on the basis of Eq. ([5]) is T(J ' ftp — > 7^ c ) = 
2.85 keV. Here we have taken the photon energy to be 114.3 MeV based on m(J/ip) = 
3096.916 MeV and m(i] c ) = 2980.4 MeV, and have assumed that the matrix element of 
jo(kr/2) between initial and final states is 1. With T to t(J/ip) = (93.4 ± 2.1) keV, this 
im plies a branching r atio B(J/ip — > 7?y c ) = (3.05 ± 0.07)%. The branching ratio observed 



in 



Gaiser et al. 



(119861 1 is considerably less, B exp (J/ijj — > 77/ c ) = (1.27 ± 0.36)%, calling for 
re-examination both of theory and experiment. 

One might be tempted to ascribe the discrepancy to relativistic corrections or the lack 
of wave function overlap generated by a relatively strong hyp erfine splitting. A calculation 



based on lattice QCD does not yet provide a definitive answer fjDudek et al 



20061 ). though it 



tends to f avor a larger decay rate . Theoretical progress may also be made using a NRQCD 



approach (jBrambilla et al 



2006). Part of the ambiguity is associated with the effective 
value of the charmed quark mass, which we take to be 1.5 GeV/c 2 . 



2. New measurements of leptonic branching ratios 



Ne w leptonic J j ip branching ratios were measured by the CLEO Collaboration (ILi et al. 
20051 ) by comparing the transitions ip{2S) — > 7r + 7r~ J/ip(lS) — > 7r + 7r~X with — 
n + n-J/^(lS) -> n+n-e+e-. The results, B(J/ip e + e~) = (5.945 ± 0.067 ± 0.042)% 
B{J/ip -> = (5.960±0.065±0.050)%, and B(J/ip -> l+f) = (5.953±0.056±0.042)% 

are all consistent with, but more precise than, previous measurements. 



3. Hadronic, gg'y, and 777 decays: Extraction of as 

The partial decay rate of J / ip to hadrons through the three-gluon final state in principle 
provides information on as{rn 2 c ) through the ratio 



T(J/ip -> ggg) 



m(J/ij)) 
2m r 



1 2 



1 + 1.6 



71 



(19) 



Both processes are governed by | \l/(0)| 2 , the sq u ared m agnitude of the 5"- wave charmonium 



wave function at the origin. In 



Kwong et all <ll988h a value of a s {m 2 c ) = 0.175 ± 0.008 



was extracted from this ratio, which at the time was measured to be 9.0 ± 1.3. This is far 
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belo w what one expects from the running of a$ down to low momentum scales (as(m 



0.3 fIBethke 



2007 



Davier et al. 



2007 



Kluth 



2006 



Kwong et al 



19881 )). highlighting the 



importance of relativistic corrections to Eq. (lT9"j) . We shall update the value of the ratio as 
extracted from data, but the qualitative conclusion will remain the same. 

The branching ratio B(J/ip — > 999) is inferred by counting all other decays, to 7%, £ + £~, 



- y* — » hadron s , and 'jgg. As mentioned earlier, we have B(J/ip — > r yr] Ci 
( Gaiser et al . 198^ ) and B(J/tp -> £+£') = (5.953 ± 0.056 ± 0.042)% ( 



1.27 ± 0.36)% 



Li et al 



20051 ) for 



£ = e,fi. We use the value R e+e - = 2.28 ± 0.04 at E cm /c 2 = m(J/V>) (ISethl . |200J) and the 
leptonic branching ratio to estimate 



B(J/ip -> 7* -> hadrons) = R e+e -B{J/i) 



;i3.6±0.3)% 



(20) 



Thus the branching ratio of J/ip to states other than ggg + gg'y is [(1.27 ±0.36) + (2 + 2.28 ± 



0.04) (5.953 ±0,070) 
(10 ±4)% (Lepage 



(26.75 ± 0.53)]%. Finally, we use V(J/if) -> ^gg)/T{J/ip -> ###) 



1983|) to infer r(J/V> -»■ ^) = (66.6 ± 2.5)%r tot ( J/^). Then 
T(J/ib^ggg) 66.6 ± 2.5 

11.2 ±0.4 



(21) 



0.188ln*nn?- Although somewhat higher than the earlier estimate, this 



r(j/v> -> £+£-) 5.953 ± 0.070 

implying a s (m c 2 ) = O.xoo_ 003 
is still far below what we will estimate from other decays, and indicates that the small 
hadronic width of the J/ip remains a problem within a nonrelativistic approach. As men- 



(Maltoni. 


2000; 


Petrelli et al.. 


1998) 



the hadronic width of the J / ip provides a useful testing ground for any approach which seeks 
to treat relativistic effects in charmonium quantitatively. The ratio 



r(J/V> -> igg) 16 a 



r(J/V> -> ffflgO 5 a s (m2) 



1 - 2.9 



a s 



TC J 



(10 ±4)% 



(22) 



itself provides infor mation on qp(m^) within a much larger range, yielding ots{m, 



0.19^q;q5 as found in lKwong et al. 



±988) 



The decay J 



sect 



i s also governed by |\P 



J/ijj — > ggg is ( iKwong e£ a/1 Il988l ) 



The ratio of its rate to that for 



T(J/tp -> 777) 54 



n \ 3 1 - 12.7a s /vr 128 ( a \ 3 1 



12.7« 5 /vr 



(23) 



T(J/ip —> ggg) 5 ^ \as J 1 — 3.7«s/7r 135 \ocs J l — 3.7a s /7i 
The last ratio is a QCD correction; eQ = 2/3 for the charmed quark's charge. (For the Y(IS') 
ratio, take eg = —1/3 and replace 3.7 by 4.9 in the denominator of the QCD correction term.) 
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With as{m 2 c ) = 0.3, the uncorrected ratio is 1.4 x 10 5 . The large negative QCD correction 
indicates that this is only a rough estimate but probably an upper bound. 



B. The r] c 



Some progress has been made in pinning down properties of the r] c (lS), but better mea- 
surements of its mass, total width, and two-photon partial width would still be welcome. 



The mass has been determined through fits t o the invariant mass spectrum o 



decay products in reaction s such as 77 — > rj c (lS) (lAsner et al 



2004 



Au 



jert et al 



2000, 



HAIS) 



2004bJ), 



200 3J) using 



B -> r] c (\S)K flFang et 0/1120031 ). and J/i/j, i/j(2S) -> 77/ c (lS) flBai et al 
all-ch arged or dominantly charged final states, and in pp — > r) c (lS) —>■ 77 (lAmbrogiani et al. 
20031 ). All these recent measurements have uncertainties in the few-MeV range, but do not 
agree wit h each other par ticularly well. The averaged value is m(rj c (lS)) = (2980.4 ± 
1.2) MeV (jYao et all 120061 ). which includes an error inflation of S — 1.5 to account for the 
spread of results. The observed splitting of 116.5 ± 1.2 MeV betwee n J/ijj and yd IS) is 



consistent with an unquenched lattice QCD prediction of ~ 110 MeV (IDavies et al 



20061 ) 



The square of th e wave function at the origin cancels out in the ratio of partial widths 
flKwong et all Il988f ). 



r(7? c -> 77) 



1 + 1.96 



Yao et al 



as_ 



J2006|), T( Vc ^ 77) 



(24) 

(7.2 ± 0.7 ± 



T(J/ip -> fi+fi- 
Using the "evaluated" partial widths in 
2.0) keV and T(J/ip -> /i + /i") = (5.55 ± 0.14 ± 0.02) keV, one finds that (3/4)r(r/ c -> 
77)/r(J/'0 — > /i + /i _ ) = 0.97±0.29, which is consistent with Eq. (|24j) but still not precisely 



enough determined to test the QCD correction. A more precise test would have taken into 
account m(J /ip) 7^ 2m c and the running of as- 

The total width of rj c is dominated by the qq fi nal st ate. Its value has not remained 



particularly stable over the years, with lYao et all (120061 ) quoting T tot (r] c ) = (25.5 ± 3.4) 
MeV. This va l ue is (3.54 ± 1.14) x 10 3 that of T(r] c 77). The gg/77 ratio is predicted 
f lKwong et all Il988l ) to be 



nvc^gg) _ 9[a 5 (m c 2 )] 2 



a z 



r(?7c -> 77) 

leading to as{m 2 c ) = 0.3w_ 05 
large QCD correction factor 1 + 8.20:5/^ ~ 1.8 



1 + 8.2^ 

7T . 



(25) 



0.30lnn 3 . This value should be regarded with caution in view of the 
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New measurements have been reported of the product of the two-photon wid ths and 
branching ratios to selected four- meson final states for the r\ c (jUehara et all 120081 ). Com- 
bining with bran ching ratios from the Particle Data Group, one obtains T(r] c — > 77) = 



2.46 ± 0.60 keV (IMetreveli 



20071 ). a value considerably lower than that just quoted, and 



disagreeing with the prediction in Eq. 
C. P-wave Xcj states 

The IP states of charmonium, Xcj, were first seen in radiative decays from the ip(2S). 
The Xcj states lie 128 / 171 / 261 MeV (J = 2/1/0) below the ij){2S). The i r masses can most 



2003 



2005b 



Armstrong et al. 



accur ately be determined in pp collisions (lAndreotti et al. 
19921 ) with Xcj jJ/ip — * 7(e + e~) or Xco — > vr°7r° by measuring the excitation curve, where 



the well-kno wn and small beam en ergy spread results in very low systematic uncertainty 
(0(100 keV (lAndreotti et all 120031 )). In principle, a precise measurement of the photon 
energy in tp(2S) — > JXcj allows a mass measurement as well, given that the ip(2S) mass 



followed by p 



loton conversions 



is very well known. BES used the decay ip(2S) 

7 — > e + e~ to improve upon the photon energy resolution (lAblikim et al . 

The J = state is wide, about 10 MeV, while the J = 1 and J = 2 states are narrower 



2005dj 



(0.89 ± 0.05 MeV and 2.06 ± 0.12 MeV, respectively flYao et all , \200m ). which is below de 



tector resolution for most exclusive Xcj decays. The most accurate width determinations 



to date come from pp experiments, again from fits to the excitation curve (lAndreotti et al. 



2003 



2005bJ) 



1. Production and decay via El transitions 

El transitions have played an important role in quarkonium physics with the initial 
theoretical papers describing charmonium suggesting that the tri plet IP states could be 



observed through the El trans i tions from the ip{2S) resonance (jAppelquist et al 



Eichten et al. 



1975 



1976 



1978 



calculations by the Cornell group (jEichten et al. 



1975 : 



rear si 
. 1975. 


iccess 
1976, 


01 xnis 
1978. 


i PICTI 

1980) 



of the present experimental values. 

New studies have been performed by the CLEO Collaboration of the rates for ip(2S) 
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TABLE III Pro perties of ik(2S) — » 'JXcJ decays, using results from lYao et al\ (j2006l ) and branching 
fractions B from lAthar et all (|2004j), as well as Eq. ([6]). 



J 




B 


lW2S)- 7XcJ ] 


\(lP\r\2S)\ 




(MeV) 


(%) 


(keV) 


(GeV- 1 ) 


2 


127.60±0.09 « 


}.33±0.14±0.61 


31.4±2.4 


2.51±0.10 


1 


171.26±0.07 « 


}.07±0.11±0.54 


30.6±2.2 


2.05±0.08 





261.35±0.33 « 


}.22±0.11±0.46 


31.1±2.0 


1.90±0.06 



7XcO,i,2 ( lAthar et all . I2004J ) and ^(25") -> 7x^,1,2 -> 11 J ( jAdam et a/1 l2005al ). We shall 
use these data to extract the magnitudes of electric dipole matrix elements and compare 
them with various predictions. 

The inclusive branching ratios and inferred rates for iJj(2S) — » lYr..j are summarized in 



Table [TTT1 Photon ene r gies ar e based on masses quoted in 



Yao et al. 



are from 



Athar et al. 



337 ± 13 keV flYao et al 



(120061 ). Branching ratios 



mm 



'artial widths are obtained from these using T tot [ip(2S)] = 
20061 ). The El matrix elements KlP^^S")! extracted using the 
nonrelativistic expression ([6]) are shown in the last column. 

In the nonrelativistic limit the dipole matrix elements in 3 S\ — > 3 Pj transitions, 
|( r )NR.|, for different J values are independent of J. Predictions of specific nonrelativis- 
tic potential models sit in a small range from 2.4 to 2.7 GeV" 1 (see Fig. [3]), with a 
slightly larger range obtained usin g potentials constructed from charmonium and bb data 



using inverse-scattering methods ( 



Quigg and Rosnerl . Il98ll ). However the magnitudes of 



the matrix elements are observed with the ordering |(x c2 |r|-0(25'))| > |(x c i|r| , 0(25')}| > 
I (■x r o\ r\i(}(2S))\. This is in accord with predictions that take into account relativistic correc 



1983 



tions ( Ebert et al 



2003a 



Moxhay and Rosner 



1985 



Grotch et al 



1984 



McClary and Byerj . 



Godfrey and Isgurl . 
1983L Figure [3] shows that at least some models are in good 



agreement with the observed rates so that we can conclude that relativistic corrections can 
explain the observed rates. However, it is probably premature to say that the transitions 
are totally understood given the large scatter of the predictions around the observed values. 

Information on the electromagnetic cascades ip(2S) — > iXcJ — > 11 J/ip is summarized 
in Table [TV] The products £>i£> 2 = B[ip(2S) — > ix c J\B\x c j ~^ iJ/^P] are taken from 
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FIG. 3 El dipole transition matrix elements for the charmonium decays 2 3 S"i — > l 3 Pj. The hor- 
izontal bands indicate the experimental results. The circles designate nonrelativistic predictions 
and the triangles relativistic predictions. Within these subsets the results are given in chrono- 



logica 



1976 



order of t he publication date. The labels refer to C-Corn ell Model (Eichtcn et al 



1978, 



1980) 



Buchmiiller Tye ( 



QR-Quigg Rosner, c c p 



Buchmuller and Tvd. 



MB-McClary Byers (jMcClarv and Bvers . 



GOS-Grotc h Owen Sebastian (IGrotch et al. 
functions of 



Godf rey and Isgur 



Galkin ( Ebert et al. 



2003aT ). 



2 and bb potentials (jQuigg and Rosnerl. 



198lh. GRR-Gupta Radford Repko (IGupta et al. 



1983), MR-Moxhay Rosner (jMoxhav and Rosner 



1981 




1984 ) , GI-Godfre y Isgur , calcu lated using the wave- 



(|1985T ). L-Lahde, DYN column (|Lahde 



2003|), EFG-Ebert Faustov 
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TABLE IV Properties of the exclusive transitions ip(2S) — > "fXcJ — > 77 •W'- 



J B X B 2 (%) B 2 (%) T tot (MeV) 

(Adam e£ a/. . 2005a ) fYao e£ a/. . , 



2 1.85±0.04±0.07 20.1±1.0 2.06±0.12 
1 3.44±0.06±0.13 35.6±1.9 0.89±0.05 
0.18±0.01±0.02 1.30±0.11 10.4±0.7 



TABLE V Properties of the transitions Xcj -> jJ/4> (|Yao et all (j20od ): Eq. ©) 



J ky 


T{ X cj^lJ/^) |(15|r|lP)| 


(MeV) 


(keV) 


(GeV)- 1 


2 429.63±0.08 


416±32 


1.91±0.07 


1 389.36±0.07 


317±25 


1.93±0.08 


303.05±0.32 


135±15 


1.84±0.10 



Adam et al. 



from 



ble (fYao et al 



(2005a ). Th ese and prior measurements may be combined with values of B\ 



Athar et al\ (120041 ) and previous references to obtain the values of B 2 in the Ta- 
2006^ . Other data com e from the high-statistics studies of Fermilab Ex- 



periment E835 (lAndreotti et all l2005bl ). who also measure total Xcj widths and present 
partial widths for Xcj — ► jJ/i>- 

The partial widths for Xcj —> lJ extracted from PDG averages for B 2 and the values of 
rtot(Xc2,i,o) mentioned above are summarized in Table PVT. The dipole matrix elements have 
be en extract e d usin g Eq. (jHJ) using photon energies obtained from the Xcj and J/ift masses 



in 



Yao et al. 



([20061). 

Predictions from both nonrelativistic and relativistic calculations are shown in Fig. |H 
Overall the nonrelativistic calculations, with typical values of 1.9 to 2.2 GeV -1 , are in 
reasonable agreement with the observed values reflecting their relative J-independence. The 
predictions including relativistic corrections are generally poorer which is surprising because 
both the IP and 15* wavefunctions have no nodes so that the integrals should be relatively 
insensitive to details of the calculation. 
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FIG. 4 El dipole transition matrix elements for the charmonium decays 1 3 -Pj — > l 3 Si. Labels are 
as in Fig. [3j 



5. Search for M2 transitions 

Attempts have been made to observe magnetic quadrupole (M2) transitions in charmo- 
nium through their interference with the dominant El amplitudes. These are not yet conclu- 
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TABLE VI Predicted and observed M2/(E1 2 + M2 2 ) 1 / 2 ratios for the transitions Xcj -> lJ/tp- 



State Prediction Experiment 

(Ambroriani et al. . 20021 fYao et al . 2006 1 

Xd -0.065(1 + k c ) -0.002t°;°°? 

Xc2 -0.096(1 + n c ) -0.13 ±0.05 



sive ( lAmbrogiani et al. 



2002 



Oreglia et all Il982l ). The best prospects a re expected for the 



most energetic photons, i.e., those in Xcj — > jJ/4>- Using the notation of ( lAmbrogiani et al. 



20021 ) . the expected normalized M2/E1 amplitude ratios a 2 for these decays are 

a 2 (Xd) = £ 71 (l + K c )/(4m c ) , (26) 
a 2 (Xc2) = (3/V5)^ 72 (l+K c )/(4m c ) , (27) 



and are shown in Tab 
in 



Oreglia et al. 



e PVTl These values are based on averages (jYao et all 120061 ) of those 



(11982( 1 and 



Ambrogiani et al. 



(J2002]). We note that a comparison between 
the ratios of the two decays would yield a more stringent test due to the cancellation of the 
charm quark mass (theory) and possible systematic uncertainties (experiment). 



3. Hadronic and 77 decays 



In principle the measured Xcj widths (lYao et al\ . 120061 ) can be used to determine 05 (m 2 ) 
if the value of the derivative of the L = 1 radial wave function for zero separation, |i?^ P (0)|, 
is known. Potential models or lattice gauge theories can be used to estimate such quantities. 
However, they cancel out in ratios of partial widths to various final states. We shall con- 
centrate on the ratios T yy ( y r j) /T gg (xr..j) for J = 2,0 (y r1 cannot d ecay into two photons). 



These are predicted to be (jEbert et al. 



2003bl ; iKwong et all Il988h 



^7l(.XcJ ) 



5cr 



-C.i ; Co 



(16a s )/(37r) 



a, 



1 + (0.2a s )/vr 



(28) 



r gg (Xcj) 9[a 5 (m 2 )] 2 ^ ' ~ z l-(2.2a s )/7r ' ^ l + (9.5a s )/n ' 
Here we have exhibited the corrections separately to the 77 partial widths (numerators) and 
gg partial widths (denominators). 

CLEO has reported a measurement of Tj xrv. ~» 77) = 559±57±48±36 eV based on 14.4 



fb- 1 of e + e~ data at ^/s = 9.46-11.30 GeV (IDobbs et alj |2006bJ ). The result is compatible 
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with other measurements when they are corrected for CLEO's B(xc2 — > jJ/ip) and B(J/ip — > 
£ + £~). The errors given are statistical, syste matic, and AB(xn 2 — > 'jJ/ip). One can average 



the CLEO measurement with a Belle result (Abe et al. 



2002a|) that is likewise corrected for 



updated input br anching fractions to obtain T(xc2 — > 77) = 565 ± 62 eV. Using T 



2.06 ± 0.12 MeV (Yao et al. 



2006J) and B{xc2 -> tW) = (20.2 ± 1.0)% flYao a/. 



tot 



[xA. 



2006|) 



one finds r(x C 2 — > gg) ~ r(x C 2 — > light hadrons) = 1.64 ±0.10 MeV. This can be compared 
to r(x C 2 — > 77), taking account of the QCD radiative corrections noted above, to obtain 

a s (ml) = 0.2961*;™ 

T he decay Xr.n ~» 77. also has been measured. Res ults from the Fermilab E835 Collabora- 



2000 



Andreotti et al. 



2004) are combined with other data to yield 
20061 ). or, with r tot (x c0 ) = 10.4 ± 0.7 MeV 



tion (lAmbrogiani et all 

B(y r() -^77) = ( 2.76 ± 0.33) x 10" 4 flYao et al. 
(IYao et all . 120061 ). r(x c0 -> 77) = (2.87 ±0 39) k eV Taking account of the (1.30 ± 0.11)% 
branching ratio of Xco t° T-^/V' ( Yao et al . 20061 ) one estimates r(x c0 — > gr^) = 10.3 ± 0.7 



MeV and hence i3(x c0 -> 77)/^(Xco -> ^) = (2.80 ± 0.42) x 10" 4 . Using Eq. (EH]) one then 
finds as{Tn 2 c ) = 0.32 ± 0.02, compatible both with the value found from the corresponding 
Xc2 rat io and with a s lightly higher value obtained by extrapolation from higher momentum 



scales ( IBethke 



2007: iDavier et al. 



20071 : 



Kluth 



20061 ) 



The success of the above picture must be regarded with some caution, as the experimental 
values of the ratios 



R 



77 



r(xc2 -> 77) 
r(xco -> 77) 



R, 



gg 



T(xc2 ->• gg) 

r(x c o -> ^) 



namely i? 77 = 0.197 ± 0.034, R gg = 0.159 ± 0.015, are far from their predicted values 



R 



1.70a, 



77 



R 



0.70a, 



15 1 + 0.06a.. 
for the nominal value as{rn 2 c ) = 0.3, which are R 



■gg 



15 1 + 3.02a a 



(29) 



(30) 



77 



0.128 and R gg = 0.110. This may be 



due to the large values of some of the first-order QCD corrections (particularly for Xc2 —>■ 77 
and Xco gg), rendering a perturbation expansion unrelia ble; it could signify effects o. 



2000 



Petrelli et al. 



negle cted color-octet components of the Xcj wavefunctions (IMaltoni 
19981 ): or it could signify that the values of |i?^ p (0)| differ for the 3 P 2 and 3 P states. It 
would be interesting to see if lattice gauge theories could shed light on this last possibility. 

The measurements of the prod uct of two-photon wid ths and branching ratios to 2(7r + 7r - ), 
K + K~7t + 7t~ , and 2(K + K~) in (jUehara et all 120081 ) for Xco,2 lead (combinin g with the 



relevant branching branching ratios from the Particle Data Group (jYao et al. 



200J)) to 
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r(x c0 -> 77) = (1.99 ± 0.24) keV and T(xc2 -> 77) = (0.44 ± 0.06) keV. The results 
entail a value i? 77 = 0.22 ± 0.04, even farther from the prediction based on first-order QCD 
corrections. 



D. The V(25) 



The ip(2S) resonance w as discovered at SLA C in e + e~ collisions within days after the 
announcement of the J / ip (jAbrams et all Il974h . 

The most precise iJj( 2 S) m ass measurement to date comes, as for the J/ip, from 
KEDR (lAulchenko et all 120031 ) . at a relative uncertainty of 7 x 10~ 6 . The current world 
average is m(ip(2S)) = 3686.093 ± 0.034 MeV. 

The total ip{2S) wi dth has been determined in direct pp production (E760 fro m the shape 



of the resonance curve ( [Armstrong et a/.l . ll993l )) as well as in e + e collisions (BES (]Bai et al. 



20021 ) from a fit to the cross-sections ip{2S) — *> hadrons, Tr + ir~J/i/), and to obtain the 

corresponding partial widths; the total width is computed as the sum of hadronic and 
leptonic widths). The PDG average of these two "direct" measurements is 277 ± 22keV. 
(Not included in the average is a recent value of 290 ± 25 ± 4 keV based on a measureme nt 
of the shape of the resonance curve by Fermilab Expe riment E835 (lAndreotti et aZ.L 120071 ) .) 
Another estimation comes from the PDG's global fit (jYao et all 120061 ) . which among many 
other measurements takes a measuremen t of r R fi into account. As for the J/ip, the radiative 
return process can be used (lAdam et all 120061 ); the decay chain presented there is e + e _ — > 
7 , 0(2S') — > j(X + J/ip), which holds for any decay ip(2S) — > XJ/ip. The observed cross- 
section is proportional to Y ee {ip{2S)) x B(J/if) — > XJ/if)), where X = 7r 7r , r] were 
used. The result of the global fit is 337 ± 13 keV. 

The two largest modern on- resonance samples are 29 M ip(2S) decays from the CLEO 
detector and a 14 M sample collected with the BES II detector. We have already discussed 
the transitions ip(2S) — > ^Xcj i n the previous subsection. Here we treat a variety of other 
electromagnetic and hadronic transitions of the ip(2S). We also briefly comment on ip{2S) 
decay via cc annihilation. 
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1. Decay to jr] c (lS) 



The decay ip{2S) — > jr) c (lS) is a forbidden magnetic dipole (Ml) transition, which would 
vanish in the limit of zero photon energy because of the orthogonality of IS and 2S wave func- 
tions. The photon energy is 638 MeV, leading to a non-zero matrix element (lSljdkr / 2)\2S) . 

I 

The decay was first observed by the Crystal Ball Collaboration (IGaiser et all 119861 ) in the 
inclusive photon spectrum of ip(2S) decays with branching ratio (2.8 ± 0.6) x 10~ 3 . The 
CLEO Collaboration measures B[ip(2S) -> yr}c{lS)] = (3.2 ± 0.4 ± 0.6) x 10" 3 , also us- 
ing the inclusive ip(2S) photon spectrum. We note that the yield fit depends consider- 
ably on the t] c width. The Crystal Ball Collaboration arrived at a width that is substan- 
tially below more recent experimental data, 11.5 ± 4.5 MeV as opposed to about 25 MeV. 
CLEO's result is for a nominal width of 24.8 ± 4.9 MeV; rescaled to the width found by 
Crystal Ball the CLEO result becomes B[ip(2S) -> 7^ C (15)] = (2.5 ± 0.6) x lO" 3 . We 
average the two primary results and arrive at (3.0 ± 0.5) x 10~ 3 . When combined with 
T tot [ip(2S)) = (337 ± 13) keV, this implies T[tp(2S) -> 7r/ c (15)] = (1.00 ± 0.16) keV, and 
hence [via Eq. ©] \(lS\j (kr/2)\2S)\ = 0.045 ± 0.004. While this result is in agree ment 



with some quark mo del predictions - for example 



Ebert et al. 



Eichten et al. 



(1975. 



1976 



1978 



19801 ) and 



2003 a|) give 0.053 and 0.042, respective 



tions ( 


Godfrev and Isgi, 


ir. 


1985; 


Grotch et al. 


1984 


Zambetakis and Bvers 


> 


1983; 


Zhang et al. 


1991 


)• 



y -, there is a wide sca t ter of predic 



Kang and Sucherl . 



1978 



Lahdd . 



2003 ; 



prediction from lattice QCD for this matrix element, as well as for corresponding forbidden 
matrix elements in the bb system. 



2. Decay to <yr] c (2S) 



The decay ip{2S) — > r yT] c (2S) is an allowed Ml transition and thus should be characterized 
by a matrix element (2S\jo{kr /2)\2S) of order unity in the limit of small k. One may estimate 
the branching ratio B[ip(2S) — > 7?7 C (2S')] b y scaling from J /ijj — > ■jtj c (1S). 

W ith B(J/^_->^r) c ) = (1.27 ±0.36)% (IGaiser et all Il986l ) and r tot (J/V>) = (93.4 ±2.1) 



keV ( Yao et al. 



20061 ). one has T{J/if) 



7?7 C ) = (1.19 ± 0.34) keV. Assuming that the 
matrix elements for ip(2S) — > 7?7 C (25') and J/tp(lS) jr} c (lS) are equal, the 2S — > 2S rate 
should be [£ 7 (2S -> 2S)/E 1 (1S -> IS)] 3 times that for IS -> IS. With photon energies of 
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TABLE VII Branching ratios for hadronic transitions ip(2S) — > J/ipX (jAdam et all l2005al ) . 

Channel B (%) 

■k+tt-J/^ 33.54±0.14±1.10 
7r°vr J/^ 16.52±0.14±0.58 
riJ/t/j 3.25±0.06±0.11 
7r°J/V> 0.13±0.01±0.01 
XJ/ip 59.50±0.15±1.90 



47.8 MeV for 2S -> 2S and 114.3 MeV for IS -> 15, this factor is 0.073, giving a predicted 



pa rtial width T 



in 



Barnes al. 



i p(2S) -> 7^ c (25)] = (87 ± 25) eV (compare, for example with 17 - 210 eV 
(booj)). Using r tot (^(25) = (337 ± 13) keV flYao et all \200& \. one then 



finds B[ip(2S) —>■ jr) c (2S)] = (2.6 ±0.7) x 10 4 , below the sensitivity of current experiments. 



3. Hadronic transitions from ip(2S) to J/ip 



The transitions ip{2S) — > n + n J/ip and ip{2S) — > -n^^J/ip are thought to proceed via 



elec tric dipole emission of a pair of gluons fo 



7T7T ( 



Bhanot et al. 



1979 



Bhanot and Peskin 



lowed by hadronization of the gluon 



1979 



Gottfried 



1978 



Peskinl . 



1979 



pair into 



Voloshinl . 



19791 ). In addition, the hadronic transitions ip(2S) — > r/J/ip and ip(2S) — > 7T J/ip have 



been observed. Rec ent CLEO measurements of the branching ratios for these transitions 



( jAdam et al. 



2005al ) are summarized in Table IVIIi (We have already quoted the branching 



ratios to J/ip via the Xcj states in Table [TVT ) 

Isospin predicts the 7r°7r° rate to be one-half t hat of 7r + 7r ~. C 



B(ti 7v J/iIj)/B(7i + 7i- J/i)) = (49.24 ± 0.47 ± 0.86)% (jAdam et all |2005a|), taking cancel- 



jEO determines 



lations of common uncertainties into account. Two other direct measurements of this 
ratio are: (57.0 ± 0.9 ± 2.6)% (BES, 



Ablikim et al. 



Andreotti et al. 



(I200M)) (57.1 ± 18 ± 4 4)% (E835 



(|2005ah ): the PDG fit result is (51.7 ± 1.8)% flYao et all 120061 ). The 7T°/r] 



ratio h as been m e asured as (4.1 ±0.4±0.1)% (CLEO. 
(BES, 



Bai et al. 



f )2004~bL . These re s ults ar e somewhat above theor e tical e xpectations, for 



Adam et al. 



(|2005aJ)) and (4.8 ±0.5)% 



example 1.6% quoted in iBai et all ( l2004bl ) based on 



Miller et al. 



19901 ). or 3.4% from 
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Ioffe and Shifmanl f 1980 



19811 ) ; iKuang et all <ll988h ; Maltmanl (J199l|) . The inclusive branch- 



ing ratio for ip(2S) — ► J/ipX, B = (59.50 ± 0.15 ± 1.90)%, is to be compared with the sum 
of known modes (58.9 ± 0.2 ± 2.0)%. Thus there is no evidence for any "missing" modes. 
The results imply B[ip(2S) — > light hadrons] = (16.9 ± 2.6)%, whose significance will be 
discussed presently. 



4- Light-hadron decays 

Decays to light hadrons proceed via annihilation of the cc pair into either three gluons or 
a virtual photon. This includes production of baryons. Such studies can receive substantial 
background due to continuum production of the same final state, e + e~ — > 7* — > hadrons. 
When interpreting the observed rate on the i/j(2S), interference effects between on-resonance 
and continuum production can complicate the picture. 

CLEO-c has collected a sample of 20.7pb _1 at ^fs = 3.67 GeV, while BES's below-^(2S') 
continuum data, 6.6 pb -1 , were taken at y/s = 3.65 GeV. At the two center-of-mass energies, 
the ip{2S) tail is of order 1/1000 [1/5000] compared to the peak cross-section for the two 
experiments (this number depends on the collider's beam energy spread). 

One expects Q = B[tp(2S) -> f]/B{J/^ -> /) to be comparable to B[^{2S) -> 
£+£-)/B(J/ip -> £+£-) = (12.4 ± 0.3)% (the "12% rule"), since light-quark decays are 
presumably governed by |\1/(0)| 2 as are leptonic decays. In fact, Q is much smaller than 12% 



(Adam et al.. 


2005b: 


Bai et al. 


2004a) 



2004aj). For example, Q(/wr)=(1.9±0.6)xl0- 



with a similar suppression for K^K^ . Many models have been brought forward to explain 
this behavior. Another interesting observation is that the Dalitz plot for the d ecay to 7r + 7r~7r° 



looks quite different for J/ip, i ) (2S), and the continuum below the i[)(2S) (jAblikim et al. 



2005bl ): In the case of the J/ip, the p bands dominate, while at the two higher energies the 
m(7T7i) distributions tend towar ds higher value s . Stud ies of ijj(2 S) — > VP states by CLEO 



(Adam et al. 



2005a) and BES (lAblikim et al. 



2004b 



2005al ) show that the 12% rule is 



much-better obeyed for VP decays forbidden by G-parity and hence proceeding via electro- 
magnetism (e.g., ip(2S) — > coir , pr), prj'). The AP (A for axial-vector) final state 6i 7r obeys 



the s caling prediction for both the charged and the neutral isospin configuration (jYao et al. 
2006k 
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Investigation of decays of the kind ip(2S) — > PP for P = 7r + , K + , and K° allow one to 
extract the relative phase and strength ratio between the ip(2S) — > qqq and iJj(2S) — > 7* 
amplitudes. This has been done by the CLEO and BES Collaborations (IDobbs et all (J2006a|) 
and references therein). 

CLEO has studied many exclusive multi-body final states of i/j(2S) (IBriere et all 120051 ) . 
several of which have not been reported before. Mode by mode, deviations from the 12% 
rule rarely amount to more than a factor of two. Moreover, the ratio of B\il>(2S) 
light hadronsl = (16.9 ± 2.6)% to B[J/i/> -> light hadrons] = (86.8 ± 0.4)% (lYao et al. 



20061 ) is (19.4 ± 3.1)%, which exceeds the aforementioned corresponding ratio for lepton 



pairs, (12.4 ± 0.3)%, by 2.3cr. The suppression of hadronic ip{2S) final states thus appears 
to be confined to certain species such as pw, K*K. 



The CLEO Col 



aboration has measured decays of i^(2S) to 



j aryon- antibaryon pairs 



( jPedlar et all 120051 ) . as has the BES Collaboration (lAblikim et al. 



2007aj). The branching 



ratios indicate that flavor SU(3) seems approximately valid for octet-baryon pair production. 
In all measured channels, the values of Q are either compatible with or greater than the 
expected 12% value. 

No clear pattern emerges, with some channels obeying the 12% rule while others fail 
drastically, and so the conclusion at this point is that the simplified picture as painted by 
the 12% rule is not adequate, and more refined models are necessary. 



E. The h^Pt] 



The h r (l l P-\) s tate of charmonium has been observed by CLEO (IRosner et al. 



2005 



Rubin et al. 



20051 ) via ijj{2 S) — > 7i°h n with h r — > 7?? r . These transitions are denoted by 



red (dark) arrows in Fig. \5\ ( 



Experiment E835 flAndreotti et al. 



Cassc 



and R 



osner 



20061 ). It has also been seen by Fermilab 
2005d ) in the reaction pp —>■ h c —+ 7^ — > 777, with 
13 candidate events. A search for the decay B ± — > h c K by the B elle Collabo r ation , 



however, has resulted only in an upper limit on the branching ratio (IFang et all 120061 ) 
B{B ± -> hJC*) < 3.8 x 10~ 5 for m(h c ) = 3527 MeV and B (h c -> 7^ c ) = 0.5. Attempts at 
previous observations are documented in lRubin et al\ (120051 ). 
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FIG. 5 Transitions among low-lying charmonium states. From lCassel and Rosnerl ((2006) 



1. Significance of h c mass measurement 

Hyperfine splittings test the spin-dependence and spatial behavior of the QQ force. 
Whereas these splittings are m(J/ip) — m(r) c ) = 116.5 ± 1.2 MeV for IS and m(ip(2S)) — 
m(j] c (2S)) = 48 ± 4 MeV for 2S levels, P-wave splittings should be less than a few MeV 
since the potential i s proportional to S 3 (r) for a Coulomb-like cc inter action. Latti c e QCD 



(IManke et al. 



2000 



Okamoto et al. 



20021 ) and relativistic potential (jEbert et all l2003al ) 
calculations confirm this expectation. One expects m(h c ) = m(l 1 Pi) ~ (m( 3 Pj)) = 
3525.36 ± 0.06 MeV. 



2. Detection in tp(2S) -> 7T°h c 



In the CLEO data, both inclusive and exclusive analyses saw a signal near (m( 3 Pj)). 
The exclusive analysis reconstructed r] c in 7 decay modes, while no rj c reconstruction was 
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FIG. 6 Left: Exclusive h c signal from CLEO (3 million ip(2S) decays). Data events correspond 



to open histogram; Monte Carlo background estimate is denoted by shaded histogram. The signal 
shape is a double Gaussian, obtained from signal Monte Carlo. The background shape is an 
ARGUS function. Right: Inclusive h c signal from CLEO (3 million ip(2S) decays). The curve 
denotes the background function based on generic Monte Ca rlo plus signa l . The dashed line shows 



the contribution of background alone. Both figures are from 



Rubin et al. 



(j2005h 



performed in the inclusive analysis. The exclusive signal is shown on the left in Fig. [61 A 
total of 19 candidates were identified, with a signal of 17.5 ± 4.5 events above background. 
The result of one of two inclusive analyses is shown on the right in Fig. [6j Combining 
exclusive and inclusive results yields m(h c ) = (3524.4 ± 0.6 ± 0.4) MeV, B1B2 = (4.0 ± 
0.8 ± 0.7) x 10~ 4 . The h c mass is (1.0 ± 0.6 ± 4) M eV below (m( 3 Pj)), at the edge of 
the (nonrelativistic) bound (jStubbe and Martini . Il99ll ) m(h c ) > (m( 3 Pj)) and indicating 
little P-w aye hyperfine splitting in ch armonium. The value of £>i£>2 agrees with theoretical 
estimates ([Godfrey and Rosnerl . I2OO2I ) of (10~ 3 ■ 0.4). 



3. Detection in the exclusive process pp — > h c 



3 7 



The Fermilab E835 Collaboration (lAndreotti et all l2005d ) studied a number of charmo 



nium resonances accessible in the direct pp channel using the carefully controlled p energy of 
the Fermilab Accumulator ring and a gas-jet fixed target. The signal of 13 events sits above 
an estimated background of 3 events and corresponds to a mass m(h c ) = 3525.8 ± 0.2 ± 0.2 
MeV. The signal strength is evaluated to be rp p B Vc7 = (10.0 ± 3.5, 12.0 ± 4.5) eV for 
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TtotOc) = (0.5, 1.0) MeV. With = 0.4 this woul d correspond to T h ^ = (25, 30) eV. 
(Kuang, Tuan, and Yan predicted Th c ^p p = 186 eV (IKuang et all Il988l ).) For comparison 
the partial widths of rj c , J/ip, Xco,i,2, and 4>(2S) to pp are roughly (33 ±11) keV, (203 ±9) eV, 



(2.25±0.25) keV 1 , (60±6) e V, (136 ±13 



ratios and total widths from 



Yao et al. 



eV, and (89 ±8) eV, where we have used branching 



(120061 ) . 



F. The r] c {2S) 



The claim by the Crystal Ball Collaboration (jEdwards et all Il982l ) for the first radial 
excitation of the r] c , the r] c (2S), at a mass of 3594 ± 5 MeV, remained unconfirmed for 20 
years. Then, the Belle Collaboration obs erved a candidate for r) c (2S) in B — > K(KsKir) 
( IChoi et all |2002j) and e + e _ — > J/ip + X (lAbe et all l2002bl ) at a significantly higher mass . 



An u pper limit on the decay i/j(2S) — > r yr] c (2S) by the CLEO Collaboration (lAthar et al. 



20041 ) failed to confirm the Crystal Ball state at 3594 MeV. The Belle res ult stimulates 



a stu dy of what other charmonium states could be produced in B decays (jEichten et al. 
2002k 



By studying its production in photon-photon collisions, CLEO fjAsner et all 12004 ) 



confirmed the presen ce of the new r] c (2S) candidate, as did the BaBar Collaboration 
( lAubert et all bo04bh . The mass of the r) c (2S) is found to be only 48 ± 4 MeV/c 2 be- 
low the corresponding spin-triplet ip(2S) state, a hyperfine splitting which is considerably 
less than the 116.5 ± 1.2 MeV/c 2 difference seen in the IS charmonium states (i.e., between 
the J/ip and the r) c (lS)). While potential models predict the ip{2S) — r] c (2S) splitting to 



be less than the J/ip — t] c splitting due to the smal 



state c ompared t o the I S state, most mode 



1994); 



Fu 



Zeng et al. 



cher 



s (e.g. 



er wavefunc t ion at t he origin for the 2S 



Ebert et al. 



fl2003ah: 



199ll ); I Gupta and Johnson! (119961 )). but not all ([Godfrey and Isgui 



Eic hten and Qu igg 



1985 



19951 ). predict a much larger splitting than what is observed. It is likely that the 



proximit y of the charm e d meson pa i r threshold, which can lowe r the ip(2S) mass by tens of 



MeV/c 2 (lEichten et al. 



2004 



2006 



Martin and Richard! Il982l ). plays an important role in 



the ip(2S) - r] c (2S) splitting. 

The width of the r] c (2S) can be determined in similar ways as that of the r] c (lS) in 



1 Using B(xcQ -> PP) = (2.16 ± 0.19) x 1CT 4 . 
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principle. In practice, the Ml photon ip(2S) — > 7mi^ c (2S') is difficult to measure well in an 
inclusive measurement rj c (2S) — > X due to background, and all channels rj c (2S) — > Y that 
can be faked by ip(2S) decay, which is several orders of magnitude more copious, will not be 
helpful in an exclusive measurement either. The only av ailable width measu rements come 
from two-photo n reactions: r = (6 3 ± 12.4 ± 4.0) Mev jAsner et all boOJ), T = (17.0 ± 



8.3± 2.5) MeV flAubert et a/ll2004bh . leading to an average of T = (14 ±7) MeV flYao et al. 



20061 ). The measurements are dominated by statistical uncertainties, which is dictated by 



the need to identify an i] c (2S) final state. To date, the only known decay modes of the 
rj c (2S) are KsK ± n T and 77. One may be led to try the same modes to which the r) c (lS) 
decays; the listed two measurements use rj c (2S) — > K s K ± it T . 

The CLEO Collaboration found that the product T(r] c (2S) -> ^)B{r] c (2S) -> K s Ktt) 
is only 0.18 ± 0.05 ± 0.02 times the corresponding product for rj c (lS). This could pose a 
problem for descriptions of charmonium if the branching ratios to KsKn are equal. More 
likely, the heavier rj c (2S) has more decay modes available to it, so its branching ratio to 
KsKn is likely to be less than that of the r) c (lS). 



G. The ^(3770) 



The t/>(3770) is primarily a 1 3 -Pi state with small admixtures of n 3 S\ states [notably 



1>(2S)] flEichten et al. 



2004 



2006 



Rosnerl . 120051 ). It is most easily produced i n e + e col- 



lision s, where it ap pears at 3770 Me V as a broad resonance (2 3 .0 ±2 .7 MeV (jYao et al. 



20061 ) ). Both Belle flAbe et all \2004 ) and BaBar flAubert et all l2006ah observed ^(3770) 
in B decay. The broadness of the state is due to the fact that decay to open charm DD 
is kinematically available and also allowed by quantum numbers. Final states involving D* 
and D s are not acces sible at this energy. T he mass and width have are most accurately 



determined in a scan. 



Ablikim et al. 



(l2006al ) achieves uncertainties of below 1 MeV for the 
mass and below 10% relative on the width. The leptonic width can be determined via a 
hadron production rate measurement as the cross-section is proportional to the coupling, 
r pP . The BES C ollaboration ha s been s tudying its decays to charmed and non-charmed final 

( l2004al )). and for t he past few years i t has been the subject 



states (see, e.g. 



Ablikim et al 



of dedicated studies by the CLEO Collaboration (IBriere et all l200ll ). 
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1. 0(3770) as a "charm factory 



The fact that 0(3770) lies so close to charm threshold [only about 40 MeV above 
2m(D )] makes it a a source of charmed particle pairs in a well-defined quantum state 
(without additional pions) in e + e~ collisions. An interesting question is whether the to- 
tal cross section <j(e + e~ — > 0(3770)) is nearly saturated by DD. If not, there could 



be noticeable non-DD decays of the 0(3770) flRosnerl . 120051 ). A CLEO measurement 
(IBesson et aJ . bo06ah . cx(0(377O)) = ( 6.38 ± O.OSjffffi) nb, appears very close to the CLEO 



value a(DD) = (6.39 ± 0.101^) nb flHe et all 120051 ). leaving litt le room for non-£>£> d e 



cays. Some question has nonetheless been raised by BES analyses (lAblikim et al. 



2006; 



I, 



2007d ) in which a substantial non-DD component could still be present. 



As a result of the difference between D° a nd D masses, the 0(3770) decays to D°D° 



He et al. 



fl2005h finds or(e + e~ -> 0(3770) 



more frequently than to D + D . For example, 
D + D-)/a(e+e- -> 0(3770) -»■ D°D°) = 0.776 ± 0.024+°;°^. This ratio reflects not only 
the effect of diff ering phase space, but also different final-state electromagnetic interactions 
( iVoloshinl . 120051 ). and is expected to vary somewhat as center-of-mass energy is varied over 
the resonance peak. 



2. Leptonic width and mixing 



The CLEO measurement of cx(0(377O)) mentioned above (IBesson et a/I l2006al ) also leads 



to a more precise value for the 0(3770) le ptonic widt 



keV. This enters into the quoted average (lYao et al. 



resul ts are (0.251 ± 0.026 ± 0.011) keV flAblikim et al. 



keV (lAblikim et al. 



1^1^(0(3770)) = (0.204 ±0.003^027) 
2006(l_of (0.242ig;^) keV. Subsequent 



2006af ) and (0.279 ± 0.011 ± 0.013) 



2007d ) from BES-II. These improvements allow a more precise estimate 



for the angle describing the mixing between ID and 25* states in 0(25') and ■0(3770): 
0(25) = -sin^llD!) + cos0|2 3 5 1 ) , 0(3770) = cos 0|l 3 £>i> + sin 0|2 3 Si) . (31) 



This mixing affects the ratio R 



^(3770)/V>(2S) 



of lepton ic widths of 0(25') and 0(37 70) and their 



predicted rates for El transitions to the \cj states ( Kuangl. 
analysis based on r ee (0(377O)) = 0.26 ± 0.04 keV ( 



200 



Rosner 



2 



Rosner 



2005|) gave 



20011 ). A previous 



12 ± 2)°, while 



the present leptonic width will give smaller errors on 0. The large present and anticipated 
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CLEO-c ^(3770) data sample will further constrain this value. A solution with negative <fi 
consistent with R^,{3,no)/ip{2S) gives an unphysically large rate for ip(2S) — > iXco- 

As noted earlier, the nonrelativistic predictions for the ip{2S) rates are generally too 
high, indicating the limitations of a nonrelativistic approach. We shall see that the pre- 
dict ed^ateJor^T 3770) — »■ IXco, which has recently been observed by the CLEO Collabora- 



tion ( Briere et al. 



20061 ). is also a factor of 2 too high in a nonrelativistic approach but is 



satisfactory when relativistic and coupled-channel effects are taken into account. 
3. V>(3770) transitions to irirJ/tp 

The rates for transitions of ^(3770) to inrJ/if) have been predicted on the assumpt ion that 
it is m ainly a .D-wave state with a small S'-wave admi xture as in the above example (jKuand . 
20061 ). (The sign convention for the mixing angle in iKuand (120061 ) is opposite to ours.) A 
wide range of partial widths, r(^(3770) — > 7r + 7r~ J/ip) = 26 to 147 keV, corresponding to 
branching ratios ranging fro m about 0. 1 % to 0.7%, is predicted. 

— > 7l + 7l^J/lp) 



The BES Collaboration (IBai et al. 
0.14 ± 



20051 ) finds #0(3770) 



(0.34 ± 



The CLEO Col 



aboration has measured a number of branching ratios for 



V>(3770) -> XJ/ij; flAdam et all 12009 ) : 23(^(3770) -> n+n-J/i;) = (0.189±0.020±0.020)%, 
25(^(3770) -> ttW/^) = (0.080 ±0.025 ±0.016)%, 25(^(3770) -> r]J/ip) = (0.087± 0.033 ± 
0.022)%, and £0(3770) -> 7r°J/V>) < 0.028%. Together these account for less than 1/2% 
of the total ^(3770) decays. In these analyses, the contribution from the tail of the if)(2S) 
decaying to the same final states has been subtracted incoherently. 



4- ^(3770) transitions to 7XcJ 



CLEO has recently reported results on ^(3770) — > iXcj partial widt h s, base d on the 
exclusive process ^(3770) — > iXci,2 11 J — * \ ll£~ *~£~ ( ICoan et all l2006bl ) and re- 



construction of exclusive Xcj decays (jBriere et all 120061 ). The results are shown in Ta- 



ble lyTTTl implying T. 7 B (^(3770) 



2005 



Eichten et al. 



IXcj) 



0(1%). Recent calculations (IBarnes et al. 



20041 ) including relativistic corrections are in good agreement with these 



measurements while nonrelativistic treatments overestimate T(ijj(3770) — > IXco)- The con- 
tribution from the tail of the ip{2S) decaying to the same final states has been subtracted 
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TABLE VIII Radiative decays ^(3770) 



Theoretical predictions of 



Eichten et al 



jXcj- energies, predicted and measured partial widths. 



(|2004l ) are (a) 



nonrela t ivisti c (c) and relativistic (d) predictions of 



without anc 



Barnes et al. 



(b) wi th coupled-channel effects; 



Rosnerl (j200lh . 



(|2005l ) : (e) shows predictions of 



Mode 


£ 7 (MeV) 


Predicted (keV) 


CLEO (keV) 




CYao et al. 2006 s ) 


(a) (b) (c) (d) (e) 


fBriere et al. 2006) 








7Xc2 


208.8 


3.2 3.9 4.9 3.3 24+4 


< 21 


7Xci 


251.4 


183 59 125 77 73 ± 9 


70 ± 17 


7Xc0 


339.5 


254 225 403 213 523 ± 12 


172 ± 30 



incoherently. 



5. "^(3770) transitions to light-hadron final states 



Several searches 



Zhu, 



2006 h Kt. Kq (lAblikim et al. 



( lAblikim et al. 



'or ^(3770) 



2007b 



Huang et al. 



(light hadrons), including VP (jAdams et al. 



2004d 



Cronin-Hennessy et all 



20oJ, and multi-body 



2006a 



20061 1 final states have been performed. No evidence 
was seen for any light-hadron ^(3770) mode above expectations from continuum production 
except for a margina lly significant branching ratio B(ip(3770) — > <prj) = (3.1 ± 0.7) x 10~ 4 
( lAdams et all l2006al ). indicating no obvious signature of noia-DD ^(3770) decays. 



H. Missing charmonium ID states 



In addition to the ^(3770), three more charmonium ID states are expected: the spin 
triplet 3 D 2 and 3 D 3 states and a spin singlet l D 2 state. All these remaining states are 
expected to be narrow. 

The masses of the remaining states are expected to be slightly above the ^(3770). Using 
the usual spin-dependent potentials we expect the ^D?} D?, 3 to lie a bout (+20, +20, 
+30) MeV, respectively, above the ^(3770) mass (IBrambilla et all [200J). The effects of 
coupling to decay channels may also produce important mass splittings. In one model these 
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additional shifts are (+37, +44, +59) MeV respectively flEichten et all 120041 . 1200O ). 

The J = 2 states ( 3 -D2 and 1 D 2 ) are forbidden by parity to decay into two pseu- 
doscalar D mesons. Hence these states are quite narrow. The principal decay modes for 
the 3 D 2 state are expected to be: radiative transitions (7 3 Pi and 7 3 p2), hadronic tran- 



sitions (t iti J / V>) , and to light hadrons (ggg). The tota 



400 keV tearnes and Godfrey! . 



2004 



Eichten et al. 



2002 



wid th is expected to be about 



20041 1 The principal decay modes 



for the spin-singlet 1 D 2 state are similar: radiative transitions (7 1 Pi), hadronic transitions 
(7nrr) r .(lS)), and ligh t hadrons (qq). Th e tota l widt h is expected to be about 460 keV 



2004 



Eichten et al. 



2002 



2004|). 



( IBarnes and Godfrey . 

Finally, the 3 D^ state has a Zweig-allowed s trong decay to DD but only in an F-wave 



( IBarnes and Godfrey! . 



2004 : 



Eichten et al. 



20041 ) . Hence the expected rate of this dominant 



decay is sma 



( IBarnes et al. 



1. For example, at a mass of 3 868 MeV this decay width is only 0.8 MeV 



2005 



Eichten et al. 



20041 . 120061 ). Thus other decay modes such as nnJ/ip 



and 7 3 P2 may be observable. 

Production rates for these remaining ID states in hadronic collisions or B meson decays 
are expected to be not significantly larger than those for ^(3770). Qualitatively, this is based 
on the assumption that the production of cc states with large relative orbital momentum is 
suppressed, and the states in question do not mix with S- or P-wave charmonium states. 



I. ^(4040), V(4160), and ^(4415) 



The '0(4040) and '0(4160) resonances appear as elevations in the measurement of R = 
o"(hadrons)/cr(/x + / u _ ). They are commonly identified with the 35* and 2D states of charmo- 
nium (Fig. ED) . Their parameters have undergone some refinement as a result of a recent 



analys is in ISethl (120051 ) . Results using initial state radiation events from Belle (lAbe et al. 



2007al ) indicate that the D*D and D*D* final states are populated throughout this en- 
ergy region, making interference effects between the resonances inevitable. BES has re- 
evaluate d earlier published data from a scan in the region 2 — 5 GeV in center-of-mass 
energy ( Ablikim et all . 2008 ) to arrive at estimates of masses, total widths, partial elec- 
tronic widths, and relative phases of ^(3770), ^(4040), ?/>(4160), and 0(4415). The analysis 
is the first to take interference between the states into account. Doing so affects especially 
the parameters extracted for the three upper states significantly. In summary, the treat- 
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ment of charmonium states above threshold from inclusive decays is not unambiguous, and 
parameters must be seen wi thin the contex t of th e method that was used to obtain them. 
Belle has reported a result (IPakhlova et all 120081 ) for the first exclusive decay of ?/>(4415), 
iff (4415) -> DD* (2460) -> DDix. Belle determines mass and total width of the ^(4415) 



from the m(DZ^2(2460)) distribution, achieving a result in agreement with ( jAblikim et al. 



2008|). 



Data taken at the ^(4040) and the ?/>(4160) can be useful to search for the 2P states 
through radiative decays ^(4160) — > 7X00,1,2(2^*) . Identifying the transition photon in the 
inclusive photon spectrum requires excellen t background suppression and is therefore a chal- 
lenge. The El branching fractions listed in (jBarnes . 2006 ) are, calculated for x c j(2P) masses 
chosen to be 2 3929/3940/3940 MeV for J = 2/1/0: 
V(4040) -> 7Xc2,i,o(2P): 0.7/0.3/0.1 x 10" 3 , 
V>(4160) -> 7Xc2,i,o(2P): 0.1/1.3/1.7 x 10~ 3 . 

The J = and J = 1 states can be distinguished since the decays Xco — > DD and 
Xd ~~ > DD* are possible but not the reverse. x C 2(2P) can decay to either, where the relative 
rate depends on the amount of phase space, which in turn depends on the mass. Exclusive 
decays to charmonium have not been observ ed, though CLEO has set upper limits on a 
number of final states involving charmonium (ICoan et all l2006al ). 



J. New charmonium-like states 



Many new charmonium states above DD threshold have recently been observed. While 
some of these states appear to be consistent with conventional cc states, others do not. Here 



we give a brief survey o 
found in 



t he new s tates a n d their possible interpretations. Reviews may be 
Godfrey! (120061 ): iRosnerl (]2006bl ): ISwansonl (120061 ). In all cases, the picture is not 
entirely clear. This situation could be remedied by a coherent search of the decay pattern 
to DD^*\ search for production in two-photon fusion and ISR, the study of radiative decays 
of ^(4160), and of course tighter uncertainties by way of improved statistical precision upon 
the current measurements. 



The motivation for this choice will become apparent in Section [IV. Jl 
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TABLE IX Summary of the X(3872) decay modes and searches. The two entries for D°D 7r° are 
both from Belle and are based on samples of 88 fb -1 and 414 fb -1 , respectively. 



final state 



X(3872) branching fraction 



reference 



ir+ir-J/tp (11.6 ± 1.9) x 10 _6 /#B+-or(3872)A-+ (> 10cJ ) (Aubert et al . 2006a) 



7r ir°J/ijj 

IXcl 



not seen 
< 0.9 x B, 



7T+7T - J /if) 



C Aubert et al . 2005c ) 
CChoi et al . 2003) 



jj/tp (3.3 ± 1.0 ± 0.3) x W- 6 /B B ^ X (3872)K+ (> 4<t) (Aubgrt oZ. . 2006c) 



(0.14 ±0.05) x £x(3872H^-j/v> (4.0a) 
r]J/iP < 7.7 x 10- 6 /^ x(387 2)^+ 

K+v-^J/i, (1.0 ± 0.4 ± 0.3) x X (3872)-> 7 r+ 7r -^ (4.3a) 

< 6 X 10- 5 /£b+^X(3872)^+ 

D+D- < 4 x W- 5 /B B+ ^ x(3872)K+ 

D°D°7T < 6 x W- 5 /B B+ ^ x(3872)K+ 

(12.2 ± 3.1±§;g) x W~ 5 /B B + 

^x(3872)_K'+ a (6.4it) fGokhroo et al . 

a Belle report the quoted number as branching fraction at the peak. They find a peak position that is 
slightly above that seen by other experiments for other X(3872) decays. 



fAhe et al . 2005a) 
f Aubert et al . 2004c) 
(AJbe et al . 2005a ) 
fAbe et al . 2004 ) 
fAbe et al . 2004) 
fAbe et al . 2004 ) 



1. X(3872) 



The X(3872), disco vered by Belle in B decays (jChoi et all 120031) and confirmed by BaBar 



(Aubert et al. 



2005d) and in hadronic production by CDF (lAcosta et all 120041 ) and DO 



( lAbazov et all 120041 ) . is a narrow state of mass 3872 MeV that was first se en decaying to 



J/'07r + 7r . No signal at this mass was seen in B — > X K, X — ► tt 7i°J/ijj ( 1 Aubert et al. 



2005d ). which would have implied a charged par tner of AT (38 7 2). It was not observed in 
two-photon production or initial state radiation (IDobbs et all 120051 ) . Subsequent studies 
focused on determining the mass, width, and decay properties in order to establish its 
quantum numbers and possible position in the charmonium system of states. To date, 
decays to -k + h~ J/ip, jJ/tp, 7r + 7r~7r° J/ip and possibly D°D 7i-° have been reported. Results 
on decay modes of X(3872) are summarized in Table HXl 



The averaged mass of this state is M = 38 71.2 ± 0.5 MeV (Yao et all 120061 ): the width 



is determined to be T < 2.3 MeV (90% C.L.) flChoi et al. 



20031 ). below detector resolution. 
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PRL91, 262001 (2003) 



PRL93, 162002 (2004) 



3900307-007 



800 D0 



X(3872) 




3.82 3.84 3.86 3.88 3.9 3.92 

M(J/\|f roc) (GeV) 



0.6 0.7 0.8 0.9 

M rV ^- - M, v (GeV/c 2 ) 



FIG. 7 Observation of X(3872) — > tt + tt J/il> in B decay (e xample from Belle (jChoi et al 



and in pp collisions (example from DO (jAbazov et al 



2004m . 



200J)) 



Signal distributions from tw o experiments are show n in Figure [TJ and mass measurements 
(including m(D°) + m(D*°) (jCawlfield et all I2007I )) are compared in Figure El 

The combined branching fraction product from Belle and BaBar is 



B[B + -> if+X(3872)] x £[X(3872) -> tt+tt- J/^] = (11.4 ± 2.0) x 10" 6 (lYao a/. 



2006) 



After s etting a limit of B[B + -> K+X(3872)] < 3.2 xlO" 4 (90% C.L.), BaBar (lAubert a/. 
2006af ) derives £[X(3872) -> tt+tt-J/V] > 4.2% (90% C.L.). For comparison, examples of 



other state s abov e open flavor threshold are £>[?/>(3770) — > tt + tt J/i/j] = (1.93 ±0.28) x 10 



-3 



(lYao et al 



2006) 



10+ 3 ( lYao et al 



(par tial width 46 keV) and limits £[t/>(4040, 4160) -> tt+tt" J/ip] of order 



20061 ) (partial widths ~ 100 keV). 



Decay into a pair of D mesons has not b een observed, and upper limits on the rate are in 



the range of a few times that for tt+tt J / if) (lAbe et al. 



2004]). A signal in B -> (D°D°tc°)K 



with m(D°D 7r ) in the right range is the first candidate for open-charm decays of X(3872). 
The observed rate is an order of magnitude above that for 7i + 7i~J/ijj. 

The dipion mass distribution favors high m(Tr+Tr~) values. This is not untypical for 
charmonium states (cf. ip{2S) — > tt+tt - J/ip), but could be an indication that the pion pair 
might even be produced in a p configuration; if that were indeed the case the X(3872) could 
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not be a charmonium state. 

The de cay X(387 2 ) — » 7r + 7r~7r° J/ij) was observed at a rate comparable to that of 
7T + 7T~J/i{) (lAbe et all . l2005al ) (preliminary). The m(7r + 7r~7r°) distribution is concentrated 
at the highest values, coinciding with the kinematic limit, which spurred speculations that 
the decay might proceed through (the low-side tail of) an u. In any case, if confirmed, 
the co-existence of both the X(3872) — > 7i + 7i^J/ip and X(3872) — > 7i + 7i^7i°J/ip transitions 
implies that the X(3872) is a mixture of both 1=0 and 1=1. 

Since the X(3872) lies well above DD threshold but is narrower than experimental res- 
olution, unnatural J p = 0~,l + ,2~ is favored. An angular distrib ution analysis by the 
Belle Collaboratio n, utilizing in part suggestions in iRosnerl (120041 ). favors J PC = 1 ++ 



flAbe et al 



2005b! ) . althou gh a higher-statistics a nalysis by CDF canno t disti n 



tween J PC = 1 ++ or 2 + (lAbulencia et all 120071 ) (see also (jKravchenkd . 



2006 



2006 



guish be- 



Marsiske 



Gokhroo et al. 



Swansonl . 120061 1 ). J PC = 2 + is disfavored by Belle's observation ( 
20061 ) of X —>■ D°D°tt°, which would require at least two units of relative orbital angular 



momentum in the three-body state, very near threshold. 

Setting aside the X(3872) — > 7t + 7t~7t°J/i/j observation for the sake of argument, among 
conventional cc states only the ID and 2P multiplets are nearby in mass. Taking into ac- 
count the angular distribution analysis, only the J PC = 1 ++ 2 3 Pi and 2~ + l 1 /^ assignments 
are possible. The decay X(3872) — > jJjjji i s observed at a rate about a quarter or less of 



2005a 



Aubert et al. 



20064 ). This would be 
and therefore 



that for X(3872) -> vr+vr-J/^ (I Abe et al . 

an El transition for 2 3 Pi but a more suppressed higher multipole for 2" 
the J PC = 1 ++ interpretation appears more likely assuming cc content. For a 1 ++ state 
the only surviving candidate is the 2 3 Pi. However, we will see that the identification of 
the Z(3931) with the 2 3 P2 implies a 2 3 P2 mass of ~ 3930 MeV, which is inconsistent with 
the 2 3 Pi inter pretation of X(3872) if the 2 3 P2 — 2 3 Pi mass splittings are decidedly lower 



than 50 MeV (IBarnes et al 



2005 



Eichten et al. 



2006h. This favors t 



re conclusion t 



X(3872) may be a P> P>°* molecule or "tetraquark" flEbert et al 



2006 



Maiani et al. 



lat the 



2005ah 



state. A prediction o f the tetraqua r 
decayin g to D D° i r (IMaiani et al. 



ration (Abe et al. 



i inter pretation is the existence of a second X particle 



2005al ). which has been reported by the Belle Collabo- 



2007al ). However, the X(3872) also has many features in common with 
tate of (D°D*° + D°D*°)/V2 ~ ccuu with J PC = 1++ Jdose and Page! . 



2004; 


Swanson . 


2004J 




iTornqvist, 


2003) 
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3900307-008 



B° J/\)/7l + 7l- + K s - 
B + -> J/\)/7l + 7C- + K + 

pp -> J/\|/7r + 7i" + X — 

PP -» J/\J/7T + 7l- + X — 
B -> J/\|/7l + 7T" + K 

average 
B^D°D°7i° + K - 
m(D°) + m(D*°) 



BaBar 

PRD 73, 011101 (2006) 

BaBar 

PRD 73, 011101 (2006) 

D0 

PRL93, 162002 (2004) 

CDF 

PRL 93, 072001 (2004) 

Belle 

PRL 91, 262001 (2003) 

PDG06 

J. Phys. G 33, 1 (2006) 

Belle 

PRL 97, 162002 (2006) 

CLEO + PDG06 

PRL 98, 092002 (2007) 
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Mass (MeV) 

FIG. 8 Comparison of mass determinations: From X(3872) — ► it + tt~ J/ip, their weighted average 
as computed by PDG, an observed threshold enhanceme nt in B — » D°D ° tt + K (at 2a deviation 



from the average), and the sum of the D° and D* mass ( Cawlfield et al. 



20071 ). 



roughly equal branching ratios is a consequence of this "m olecular" assignment. A new 
measurement of m(D°) = 1864.847 ± 0.150 ± 0.095 MeV/c 2 Jcawlfield et all 120071 ) implies 
m(D°Z>o) = 3871.81 ± 0.36 MeV/c 2 and hence a binding energy of 0.6 ± 0.6 MeV (see 
also Fig. [8]). Irrespective of its eventual interpretation, the evidence is mounting that the 
X(3872) is not a conventional cc state. 



2. Z(3930) 



Belle has reported a candidate for a 2 3 P2(Xc2(2P)) state in 77 collisions (lUehara et al. 



20061 ). decaying to DD. The state appears as an enhancement in the m(DD) distribution 
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FIG. 9 Belle's Xc2^P) candidate (|Uehara et al 



2006I ): Left: The invariant mass m{DD) distri- 
bution in two-photon production of the ^(3930), D + D~ and D°D° combined. The signal yield is 
64 ± 18 events. The two curves are fits with and without a resonance component. Right: cos#*, 
the angle of the D meson relative to the beam axis in the 77 center-of-mass frame for events with 
3.91 < m(DD) < 3.95 GeV; the data (circles) are compared with predictions for J = 2 (solid) and 
J = (dashed). The background level can be judged from the solid histogram or the interpolated 
smooth dotted curve. 



at a statistical significance of 5.3cx. The relative D + D~ and D°D° rates are consistent with 
expectations based on isospin invariance and the D + — D° mass difference. Combining 
charged and neutral modes, a fit shown in the left-hand panel of Fig. [9] yields mass and 
width M = 3929 ± 5 ± 2 MeV and T = 29 ± 10 ± 2 MeV. Although in principle the D-pair 
could be produced from D*D, the observed transverse momentum spectrum of the DD pair 
is consistent with no contribution from D*D. 

The observation of decay to DD makes it impossible for Z(3930) to be the rj c (3S) state. 
Both x c0 (2P) and x c2 (2P) are expected to decay to DD (x cl (2P) is not; it only decays to 
D*D). To distinguish between the two remaining hypotheses, the distribution in 8*, which 
is the angle of the P-meson relative to the beam axis in the 77 center-of-mass frame, is 
examined. This distribution is consistent with sin 4 6* as expected for a state with J = 
2, A = ±2 (right-hand panel of Fig. [9]). The two-photon width is, under the assumption of 
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a tensor state, measured to be T 77 • B D fj = 0.18 ± .05 ± 0.03 keV . 



BaBar has searched for Z(3930) decay into "/JM (jAubert et a/.U2006d j. and set an upper 
limit B(B -> Z(3930) + K) x £(Z(3930) -> 7J/VO < 2.5 x lO" 6 . 

The predicted mass of the x c2 (2P) is 3972 MeV and the predicted partial widths and 



2006f ) ; 



total width assuming M[2 3 P 2 (cc)] = 3930 MeV are (lEichten et al. 

r(x c2 (2P) -> £>£>) = 21.5 MeV, 

r(x c2 (2P) -> DD*) = 7.1 MeV and 

r t otai(x c2 (2P)) = 28.6 MeV, 
in good agreem ent with the e xperimental measurement. Furthermore, using r(x c2 (2P) 



77) = 0.67 keV (IBarnesl . ll992[ ) times £(x c2 (2P) -> DD) = 70% implies T iY B dD = 0.47 keV, 



which is within a factor of 2 of the observed number, fairly good agreement considering the 
typical reliability of two-photon partial width predictions. 

The observed Z(3930) properties are consistent with those predicted for the Xc 2 (2P) 
2 3 P 2 (cc) state. So far, the only mild surprise is the observed mass, which is 40 — 50 MeV 
below expectations. Adjus ting that, all other properties observ ed so far can be accomodated 



within the framework of (lEichten et al. 



2006 



Swansonl . 120061 ). The Xc 2 (2P) interpretation 



could be confirmed by observation of the DD* final state. We also note that the x c2 (2P) 
i s predicted to undergo radiative transit ions to ip{2S) with a partial width of (9(100 keV) 



( IBarnes et al. 



2005 



Eichten et al. 



3. ^(3940) 



The V(394 0) was first seen 



7r + 7r"7r°)J/y> (I Abe et al. 



yy Belle in the uoJ/%p subsystem in the decay B — > Ku(- 



2005a ). The final state is selected by kinematic constraints that 



incorporate the parent particle mass m(B) and the fact that the P-meson pair is produced 
with no additional particles. Background from decays such as i£i(1270) — > ujK is reduced by 
requiring m(u)J/ip) > 1.6 GeV. The KujJ /ip final state yield is then further examined in bins 
of m{ujJ/ip). A threshold enhancement is observed, shown in Figure [TUl which is fit with a 
threshold function suitable for phase-space production of this final state and an ^-wave Breit- 
Wigner shape. The reported mass and width of the enhancement are M = 3943±11±13 MeV 



3 Barnes, Godfrey and Swans on ( Barnes et all I2005T ) obtain similar results when the 2 3 P2 mass is rescaled 
to 3930 MeV. See lSwansor] (|2006t ). 
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FIG. 10 Belle's X c i(2P) candidate (jAbe et al 



2005d ): The invariant mass m(ujJ /ip) distribution 



in m(B — > KuiJ/t/j) decay. The signal yield is 58 ± 11 events. The two curves are fits without (left) 
and including (right) a resonance component. 



and T = 87±22±26 MeV. A fit without a resonance c ontribution gives no g ood description of 
the data. BaBar confirmed the existence of the state (lAubert et all l2007bl ). also for charged 



and neutral P decays; the values are M = 3914.6l3;4 ± 1.9 MeV and T = 33Z 8 ± 5 MeV, 
somewhat different from Belle's. 

The mass and width of V(3940) suggest a radially excited P-wave charmonium state. 
The combined branching ratio is B(B -> KY) ■ B(Y -> uJ/ip) = (7.1 ± 1.3 ± 3.1) x 1(T 5 . 
One expects that 8{B -> Kx c j{2P)) < B(B -> K X cj) = 4 x 10" 4 . This implies that 
B(Y — > uJ/ip) > 12%, which is unusual for a cc state above open charm threshold. 

For the y r i (2P ) we ex pect DD* to be the dominant decay mode with a predicted width of 



140 MeV (iBarnesl . 



20061 ). which is consistent with that of the V(3940) within the theoretical 



and experimental uncertainties. Furthermore, the \ci is also seen in P-decays. Although 



the decay 



u)J is unus ual, the corresponding decay xm(2P) — > toT(lS) has also 



been seen ( iSeverini et all , 120041 ) . One possible explanation for this unusual decay mode is 
that rescattering through DD* is responsible: 1 ++ — > DD* — > uJ/ip. Another contributing 
factor might be mixing with the possible molecular state tentatively identified with the 
V(3872). 
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BaBar has searched for F(3940) decay into 'yJ/ip (lAubert et all l2006d ). and set an upper 
limit B[B -> V(3940) + K) x £(V(3940) -> jJ/ip) < 1.4 x 10~ 5 . 

The Xci(2-P) assignment can be tested by searching for the DD and DD* final states 
and by studying their angular distributions. With the present experimental data, a Xco(2-P) 
assignment cannot be ruled out. 

4. Charmonium in e+e" -> + X: X(3940) and X(4160) 

Belle studied d ouble-charmonium production and e + e~ — > J/ip + X near the T(4S) 
( I Abe a/.l . l2007bl ) and observed enhancements for the well-known charmonium states tj c , 
Xco, and rj c (2S), at rates and masses consistent with other determinations. In addition, 
a peak at a higher energy was found. The mass and width were measured to be M = 
3936 ± 14 ± 6 MeV and T = 39 ± 26 (stat) MeV. 

To further examine the properties of this enhancement, Belle searched for exclusive decays 
J/ip — ► DD(*\ given that these decays are kinematically accessible. The J/ip recoil mass for 
the cases DD and DD* are also shown in Figure [TTl An enhancement at the X(3940) mass 
is seen for DD*, but not for DD. The mass and width determined in this study are M = 
(3943 ± 6 ± 6) MeV, T < 52 MeV (90% c.L). Note that the inclusive and exclusive samples 
have some overlap, and thus the two mass measurements are not statistically independent. 
The overlap has been eliminated for the branching fraction determination. A signal of 5.0a 
significance was seen for DD*, but none for DD. In addition, the V(3940) did not show a 
signal for a decay uJ/ip, unlike the V(3940). These findings are summarized in Table IXl 



Belle updated their study with slightly higher luminosity (I Abe et all l2007d ). The study 
confirmed the observation in the exclusive decay with comparable parameters of the V(3940) 
but higher significance, and added the following pieces of information: (1) There is no indica- 
tion of an V(3940) signal in the invariant mass spectrum of DD, but there is a statistically 
significant population spread out over a wide range. It is mandatory to understand this 
before it is possible to quantify an upper limit on X(3940) — > DD. (2) In the final state 
D*D* , a peak of 5.1a statistical significance is fit with a Breit-Wigner shape and claimed as 

a resonance of mass M = (4156^20 ± 15 ) MeV > and width r = (^-li 1 ± 21) MeV, distinct 
from the V(3940) (preliminary). 

If confirmed, the decay to DD* but not DD suggests the V(3940) has unnatural parity. 
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FIG. 11 Belle's A(3940) (|Abe et all . l2007bh . sighted in e+e" -> + A: Left: The mass of 
the system recoiling against the J/i/j. The excess at A(3940) contains 266 ± 63 events and has 
a statistical significance of 5.0a. Right: Study of A(3940) decay into D mesons, e + e~ — > J ftp + 
DD(*\ Top: DD, no signal is seen at 3940 MeV. Bottom: DD*, the signal amounts to 24.5 ± 6.9 
events (5. Oct). 



TABLE X Properties of the A (3940) (jAbe et al 



2007bh 



Mass 

Total width 
£(A(3940) -> DD*) 

£(A(3940) -> DD) 
i3(A(3940) -» wJ/V>) 



3936 ± 14 ±6 MeV (incl.) 
3943±6±6MeV (DD*) 
< 52 MeV 
(96±|±22)%, 
> 45% (90% C.L.) 

< 41% (90% CL) 

< 26% (90% CL) 



The lower mass rj c (lS) and rj c (2S) are also produced in double charm production. One 
is therefore led to try an T}J3S} ass i gnme nt, although this state is expected to have a 



somewhat higher mass flBarnes et al 



200 5J). The pre dicted width for a 3 1 S state with a 



mass of 3943 MeV is ~ 50 MeV (lEichten et all 120061 ) . which is in not too bad agreement 
with the measured AT (3940) width. 

Another possibility due to the dominant DD* final states is that the X(3940) is the 
2 3 Pi(cc) Xi(2P) state. It is natural to consider the 2P(cc) since the 2 3 Pj states are predicted 
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to lie in the 3920-3980 MeV mass region a nd the widths are predicted to be in the range 



r(2 3 Pj) = 30-165 MeV (IBarnes et a/1 120051 ) . The dominant DD* mode would then suggest 



that the X(3940) is the 2 3 Pi(cc) state. The problems with this interpretation are (1) there 
is no evidence for the l 3 Pi(cc) state in the same dat a, (2) the pred icted width of the 2 3 Pi(cc) 
is 140 MeV (assuming m(2 3 Pi(cc)) = 3943 MeV) (IBarnesl . 120061 ) . and (3) there is another 
candidate for the 2 3 P 1 (cc) state, the V(3940). 

The most likely interpretation of the V(3940) is that it is the 3 1 Sq(cc) r) c (3S) state. Tests 
of this assignment are to study the angular distribution of the DD* final state and to observe 
it in 77 — > DD*. 



5. it + it-J/4; in Initial State Radiation: Y(4260) and X(4008) 



Perhaps the most intriguing of the recently discovered states is the V(4260) reported by 
BaBar as an enhancement in the ttttJ /ip subsystem in the radiative return reaction e + e _ — > 



7isr J 1 1/J7171 ( jAubert et all l2005bl ). where "ISR" st ands for " i nitia 



state radiation." This 



and subsequent independent confirmation signals (jHe et all 120061 ; lYuan et all 120071 ) are 
shown in Fig. [121 The measured mass, width, and leptonic width times B(Y — > J/i/)7r + 7T~) 
are summarize d in the first row of Table IXII Further evidence was seen by BaBar in B — > 
K(tt + tt-J/iP) (jAubert et all \20064 ). 



The CLEO Collaboration h as confirmed the V(4260), both in a direct scan (ICoan et al. 



2006al ) and in radiative return (jHe et all 120061 ). Results from the scan are shown in Fig. 
including cross-section increases at E cm = 4260 MeV consistent with V(4260) — > 7r + 7r~ J/ip 
(llcr), tt°ti J/iI) (5.1a), and K + K~ J/ip (3.7a). There are also weak signals for ■0(4160) — > 
TT+Tv- J/ip (3.Qa) and 7r 7r°JM (2 .6a), consistent with the V(4260) tail, and for 0(4040) -> 



7T+7T- J/ip (3.3a). 
row of Table IXII 



He et al. 



(]2006l ) determines the resonance parameters shown in the second 



Belle ( jYuan et all 120071 ). also in ISR events, fit the 7r + 7r~J/0 enhancement with an 
additional component, two coherent Breit-Wigner functions in total, in order to achieve a 
better description of the low-side tail of the V(4260). The fit results in mass and width of 
M = (4008±40l2g) MeV and T = (226±44±^) MeV for the lower resonance. The values for 
the upper (the V(4260)) are listed in Table IXll Interference leads to a two-fold ambiguity in 
the rate, corresponding to constructive and destructive interference. Both solutions arrive 
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FIG 12 Y(4260) signal in ISR from the T(45) by BaBar (lAubert et al 



20061 ). and Belle (|Yuan al. 



2005b|), CLEO (|He al, 



20071 ) . The fit parameters are given in Table IXII 



at the same fit function. The two solutions differ markedly. The lower-lying state is not 
associated with any presently known charmonium state. 

The invariant mass distribution m(7r + 7r~) looks quite different for events at ~ 4260 MeV 
than above a nd below; the dist ribution is shifted towards higher values, not consistent with 

2007L 



phase space (lYuan et al. 



A variety of ratios between c hannels have b een meas ured now (lAubert et all l2006d ; 



Coan et all. 



2006a 



Gowdyl . 



2006 



Heltsley 



2006 



Ye 



20061 ). which should help narrow down 



the possible explanations of Y(4260). They are listed in Table IXlTl The preliminary upper 
limit for the ratio of DD to 7i + 7i^J/ip of 7.6 may not seem particularly tight at first glance, 
but is to be compared, for example, with the same ratio for the ^(3770), where it is about 
500. 
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TABL E XI Co mparison of pa rameters of y(4260) as mea sured by the BaBar (jAubert et all 
2005bJ), CLEO (|He et oflbood ). and Belle (|Yuan et all boOTT ) Collaborations. 



Collab. 


Mass 


r 


T ee x £(Y(4260) -» tt+ti- J/ip) 




(MeV/c 2 ) 


(MeV/c 2 ) 


(eV) 


BaBar 


4259 ± 8±g 


88 ± 23±2 


5.5 ± 1.0±g;f 


CLEO 


4284+^ ± 4 


73l| ± 5 


8.9±li ± 1.8 


Belle 


Two-resonance fit: 






4247 ± 121^ 


108 ± 19 ± 10 


6.0 ± 1.2±£| or 20.6 ± 2.3±?;^ 




Single- 


resonance fit: 






4263 ± 6 


126 ±18 


9.1 ± 1.1 



A number of explanations have app eared in the lite rature: it) {AS) ( 



cscs t etraquark (IMaiani et al. 



2005 



Zhul. 



2005bJ) , and cc hybrid (I Close and Pagd . 12005 



Janes- Estrada, 



2005), 



( iLlanes- Estrada 



_, iKou and Pend . 

200^). In so me models the mass of the F(4260) is consistent with the 4S(cc) level 
20051 ) . Indeed, a 4S ch a rmon ium level at 4260 MeV/c 2 was anticipated 



on exactly this basis (IQuigg and Rosnerl . 



19771 ). With this assignment, the nS levels of 



charmonium and bottomonium are remarkably congruent to one another. However, other 
calculation s using a linea r plus C oulomb potential identify the 4 3 S*i(cc) level with the ■0(4415) 



state (e.g. 



Barnes et al 



( 120051 )). If this is the case the first unaccounted-for 1 (cc) state 
is the ip(3 3 Di). Quark models estimate its mass to be m(3 3 Di) ~ 4500 MeV which is much 
too heavy to be the F(4260). The F(4260) therefore represents an overpopulation of the 
expected 1 states. The absence of open charm production also argues against it being a 
conventional cc state. 

The hybrid interpretation of F(4260) is appealing. The flux tube model predicts that 



the lowest cc hybrid mass is ~ 4200 MeV ( 



Barnes et al. 



having similar expectations ( jLacock et all 119971 ) . Models of hybrids typically expect the 



19951 ) with lattice gauge theory 



wavefunction at the origin to vanish implying a small e + e width in agreement with the 
observed value. Lattic e gauge theory foun d that the bb hybrids have large couplings to 



closed flavor channels (IMcNeile et all |2002j) which is similar to the BaBar observation of 
Y — > J/iIj7t + tt~; the branching ratio of B(Y — > J/ilm+n^) > 8.8% combined with the 



54 



>10h 
o 

2 




CO 



c 

> 

LU 



I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

V 




■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
• . n + n J A|/ 

m±4* ■ ■ ■ i ■ ■ 1 1 


i i i i i i i i i 


■ i i i i i 

f 1 ^ 

■ 1 htV. I 


7r 7i° J/\|/ 
' i VTttI i i i i i i i i 


1 ■ ■ ■ H ■ ■ ■ 
i . . . . i . . . 


I 

: ; 


I K + K'JA|/ . 

i; 

i — ■ i ■ ■ ■ ■ i ■ ■ ■ ■ i 



:0.15 -0.1 -0.05 0.05 0.1 0.15 

Missing Momentum (GeV) 



FIG. 13 Evidence for y(4260) from a direct scan by CLEO (jCoan et al 



2006ah 



observed width implies that T(Y — > J/?/>7r + 7r~) > 7.7±2.1 MeV. This is much larger than the 
typical charmonium transition widths of, for example, r(-0(377O) — > J/ip7r + 7T~) ~ 80 keV. 
And the Y is seen in this mode while the conventional states ^(4040), '0(4160), and ■0(4415) 
are not. 

One predicted consequence of the hybrid hypothesis is that the dominant hybrid char- 
monium open-charm decay modes are expected to b e a meson pair with an S- wave (D, D*, 
D s , D*) and a P-wave (Dj, D s j) in the final state (IClose and Pagd . 120051 ). The dominant 
decay mode is expected to be DD\ + c.c. (Subsequently we shall omit "+c.c." in cases 
where it is to be understood.) Evidence for a large DD\ signal would be strong evidence 
for the hybrid interpretation. A complication is that D pi thresho l d is 4 287 MeV/c 2 if we 
consider the lightest D\ to be the narrow state noted in lYao et al\ (120061 ) at 2422 MeV/c 2 . 
The possibility also exists that the F(4260) could be a DD\ bound state. It would decay to 
DttD*, where the D and 7r are not in a D*. Note that the dip in R e + e - occurs just below 



DDi threshold, which may be the fir st S- wave meson pair accessible in cc fragmentation 



(IClose and Pagd . 



2005 



Rosnerj, 



2006al ). In addition to the hybrid decay modes given above, 
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TABLE XII Experimental results on Y(4260) decay. The last co l umn g ives the relative rat e 



compared 
except (a) 



to 7r + 7r J/ij) for each cha nnel. Data ar e from 



Aubert et al. 



(|2006dh . (b) 



Coan et al 



Gowdj (1200J ), and (c) 



Ye ( 



(2006a 


) and 


Heltslev 


(2006) 



2006). Unless indicated otherwise, 



upper limits are at 90% CL. 



Channel 


cross-section (pb) 




7T + 7T~ J/tJj 


± 4 


1 




51 ± 12 (a) 


1 


TT°TT°J/%b 


23±§ 2 ± 1 


0.39ig;^ ± 0.02 


K+K-J/ip 


911 ±1 


0.15ig;j$ ±0.02 


7]J/lp 


< 32 


< 0.6 


TT°J/4> 


< 32 


< 0.2 


ri'J/iP 


< 19 


< 0.3 


1T + 1T~ IT® J / Ip 


< 7 


< 0.1 


7]7] J /ip 


< 44 


< 0.8 


7r+7r"V(25') 


< 20 


< 0.3 


#(2S) 


< 25 


< 0.4 


^XcO 


< 234 


< 4.0 


IXcl 


< 30 


< 0.5 


!Xc2 


< 90 


< 1.6 




< 46 


< 0.8 


TT + TT'TT°Xc2 


< 96 


< 1.7 




< 5 


< 0.1 (also see (b)) 


DD 




< 7.6 (95%CL) (c) 


PP 




< 0.13 (a) 



lattice gauge theory suggests that we search for other closed charm modes with J PC = 1 : 
J/ipr), J/iprj', Xcj^ and more. Distinguishing among the interpretations of the Y(4260) will 
likely require careful measurement of several decay modes. 

If the F(4260) is a hybrid it is expected to be a member of a multiplet consisting of eight 
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states with masses in the 4. to 4.5 GeV ma s s ran ge with lattice gauge theory preferring 



the higher end of the range (iLiao and Mankei . 12002 ). It would be most convincing if some 
of these partners were found, especially the J PC exotics. In the flux-tube model the exotic 
states have J PC = +_ , and 2 +_ while the non-exotic low-lying hybrids have 0~ + , 
2-+, 1++, and 1". 



6. States decaying to tt + tt ip(2S) 



In the radiative return process e + e — > 7 + X, the BaBar Collaboration (lAubert et al. 



2007al ) reports a broad structure decaying to 7r + ir~ip(2S), where ip(2S) — > 7i + 7i~J/tp. A 
single-resonance hypothesis with m(X) = (4324 ± 24 (stat)) MeV and T(X) = (172 ± 
33 (stat)) MeV is adequate to fit the observed mass spectrum. 

Belle, with more than twice the sample size used in the BaBar analysis, observed two 
enhancements in the same reaction (IWang et all 120071 ): One that confirms BaBar 's mea- 
surement, at M = (4361 ± 9 ± 9) MeV with a width T = (74 ± 15 ± 10) MeV (statistical 
signficance 8a), and a second, M = (4664 ± 11 ± 5) MeV with a width of T = (48 ± 15 ± 3) 
MeV (5.8cx). The existence of the higher-energy peak is not excluded in the BaBar data. 

Given the uncertainty in the masses and widths of the lower state decaying into 
7r + 7r~ip(2S) and the V(4260), the possibility that they are different manifestations of the 
same state cannot be excluded. 



V. BOTTOMONIUM 
A. Overview 

Some properties and decays of the T (bb) levels are summarized in Fig. [2j The measured 
masses of the T states below open flavor threshold have accuracies comparable to those in 
charmonium since similar techniques are used. Experimentally, the situation is more difficult 
due to the larger multiplicities involved and due to the increased continuum background 
compared to the charmonium region. 

Modern data samples are CLEO's 22M, 9M, 6M T(l,2,35') decays (with smaller off- 
resonance samples in addition) and Belle's T(3S') sample of 11M T(3S'). 
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FIG. 14 From lArtuso et al 



(|2005l ): Inclusive photon spectrum in T(nS) — ► jX, for n = 2 (left) 
and n = 3 (right), before and after background subtraction. In the upper plots, the dashed 
line indicates the background level; in the lower plots, the fit contribution for each resonance is 
delineated. The low-lying solid curves in the lower right plot show two background contributions. 
The three peaks corresponding to the x&j(l, 2P) are clearly visible. The peak position determines 
the Xfej(l>2P) masses. The signal area is used to determine the Y(nS) — ► 7x&j(l,2P) branching 
fraction. 



The x&j(l ; 2P) states are reached through El transitions; branching fractions for n — > 
n — 1 range from 4-14%. Their masses are determined from the transition photon energies. 
Their intrinsic widths are not known. E xamples of fits to th e inclusive photon spectrum that 
led to Xbj(lj2P) mass determinations (lArtuso et all 120051 ) are shown in Fig. [TH Exclusive 
hadronic decays of the Xw(l> 2P) states have not been reported; information exists only on 
transitions within the bottomonium spectrum. An T(1D) candidate has been observed; the 
singlets %(1, 2S) and hb(l, 2P) have thus far escaped detection. 

Mass differences withi n the bottomo nium spectrum are in agreement with unquenched 



lattice QCD calculations (ILepage 



20051 ). Direct photons have been observed in IS, 2S, and 



3S decays, leading to estimates of the strong fine-structure constant as consi stent with others 



Besson et al 



2006bl ). The transitions x&(2P) — > Tinxbi^P) have bee n seen ( 



Cawlfield et al. 



20061 ). BaBar has observed T(4S) -> n + n-T(lS,2S) transitions flAubert et all . l2006bh . 
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while Belle has seen T(45) -> 7r+7r-T(lS) dSokolov et o/.Ll2007f ). 

Decays to light hadrons proceed, as in the case of the charmonium states, via annihilation 
of the heavy quarks into ggg, gg^f or 7*, which subsequently hadronize. At higher energies, 
fragmentation into low-multiplicity states is suppressed, and so the second step makes it 
difficult to arrive at a simple scaling prediction to translate bottomonium and charmonium 
results into each other. Comparing the T states with each other, for example by constructing 
a prescription akin to the 12% rule in charmonium, is possible, but to date only a few 
exclusive radiative decays to light mesons, but no exclusive non-radiative decays to light 
mesons, have been observed. 



B. T(15,25,35) 



1. Masses and total widths 



The best measurements of the narrow T(nS) states, as was the case for the J PC = 1 
states in charmonium, come from fits to the cross-section T(nS) —+ hadrons around the res- 
onance together with a very precise beam energy calibration using resonant depolarization. 
This leads to precision mass determinations with uncert ainties of order 1 00 keV . 

The Tflff) mass me asurements from CUSB (IMacKay et all 11984 ) and MD- 



1 (jArtamonov et al. 



20001 ) each have a relative precis ion of one part in 10 5 , b ut are about 



0.5 MeV apa rt. The T (2S ) deter minations by MD-1 (jArtamonov et al. 



2000) and DORIS 



experiments (IBarber et all 11984) agree well. There is only one measurement of m[T(3S)), 
again by MD-1 (jArtamonov et all 120001 ) . 

The below-flavor T(nS) states are narrow, some ten keV, whereas the T(4S), for which 
the decay to BB is kinematically possible, has a full width three orders of magnitude higher. 
The intrinsic widths of the Y(l, 2,3.5) cannot be determined directly in e + e~ collisions as 
they lie well below the beam energy spread. They can be determined indirectly, by using 
the relation 

r„„ r 

(32) 



r = — 



B 



1111 



where the last step assumes lepton universality. Expanding by the hadronic partial width, 
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Thad = (1 — 3B^)/Y, the equation reads: 

T-, r ee r had /r 



(33) 



The expression in the numerator is directly accessible in the reaction e + e~ — > T(nS) — > 
hadrons; the integral of the hadronic cross-section over the resonance is proportional to the 
product of widths. The muonic branching fraction can be determined from a measurement 
of £ = r wt /Thadrons, which is independent of the total width; B m = T^/T = r A1At /(r had + 
3r w ) = £/(l + 3£). The current status of experimental precision is below 2% for r ee rh a d/r 
and 3-4% for B m . The corresponding measurements are discussed in Section IV. B. 21 



2. Leptonic branching ratios and partial widths 



New values of B\T(1S, 2S,3S) ->■ u 



= (2.49 ± 0.02 ± 0.07, 2.03 ± 0.03 ± 
0.08, 2.39 ± 0.07 ± 0.10)% (lAdams et all 120051 ). when combined with new measurements 
T„.(1S, 2S, 3S ) = (1 .354 ± 0.004 ± 0.020, 0.619 ± 0.004 ± 0.010, 0.446 ± 0.004 ± 0.007) keV 



flRosner et al. 



2006|), imply total widths T tot (lS, 2S,3S) = (54.4±1. 8, 30.5±1.4, 18.6±1.0) 



keV. The values of r tot (2S', 3S) changed considerably with respect to previous worl d av 



erages. Combining with previous data, the Particle Data Group (lYao et al. 



2006 now 



quotes T tot (lS,2S,3S) = (54.02 ± 1.25, 31.98 ± 2.63, 20.32 ± 1.85) keV, which we shall 
use in what follows. This will lead to changes in comparisons of predicted and observec 
trans ition rates. As one example, the study of T(2S,3S) — > decays (jArtuso et al. 



20051 ) has provided new branching ratios for El transitions to Xbj(^P), Xbj(2P) states. 
These may be combined with the new total widths to obtain updated partial decay widths 
[line ( a) in Table IXIII | . which may be compared with one set of nonrelativistic predic- 
tions (IKwong and Rosnerl . Il988l ) [line (b)]. The suppression of transitions to J = states 



by 10-20% with respect to nonrelativistic expectations agrees with re 



1983 



Moxhay and Rosner 



1983 



Skwarnicki 



ativistic predictions 



20051 ). The partial width 



(IMcClary and Byersl . 

for T(3S*) — » 7l 3 -Pn is found to be 61 ± 23 eV, about nine times the highly-suppressed value 



predicted in IKwong and Rosnerl (119881 ) . That prediction is very sensitive to details of wave 
functions; the discrepancy indicates the importance of relativistic distortions. 

The bra nching ratios B\T(1S , 2S, 3S) — > t + t~] have been measured by the CLEO Col- 



laboration (IBesson et al. 



20071 ). and are shown in Table IXIV1 They are consistent with 
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TABLE XIII Comparison of observed (a) and predicted (b) partial widths for 2S — ► lPj and 
35" — > 2Pj transitions in bb systems. 





r (keV), 25" — » lPj transitions 
J=0 J=l J = 2 


r (keV), 35 — * 2Pj transitions 
J=0 J=l J=2 


(a) 
(b) 


1.20±0.18 2.22±0.23 2.32±0.23 
1.39 2.18 2.14 


1.38±0.19 2.95±0.30 3.21±0.33 
1.65 2.52 2.78 



TABLE XIV Ratio R TT = B[T(nS) -» rr]/jB[T(n5) -> H and B[T(nS) -> rr] dBesson e* aZ.I . 
20071 ). 



£[T(nS) 



rr 



T(15) 1.02 ± 0.02 ± 0.05 2.54 ± 0.04 ± 0.12 
T(25) 1.04 ± 0.04 ± 0.05 2.11 ± 0.07 ± 0.13 
T(35) 1.05 ± 0.08 ± 0.05 2.52 ± 0.19 ± 0.15 



lepton universality and represent the first measurement of the T(35) — > tt branching ratio. 

3- igg/ggg ratios 

Th e direct phot o n spec trum in IS, 2S, 3S decays has been measured using CLEO III 
data (IBesson et all l2006bl ). The ratios P 7 = B(gg^) /B(ggg) are found to be Ry(lS) = 
(2.70±0.01±0.13±0.24)%, P 7 (2S) = (3.18±0.04±0.22±0.41)%, P 7 (3S) = (2.72±0.06± 
0.32 ± 0.37)%. P 7 (1S) is consistent with an earlier CLEO value of (2.54 ± 0.18 ± 0.14)%. 



C. El transitions between Xbj{fiP) and S states 

We have already discussed the inclusive branching ratios for the transitions T(2S) — > 
7Xbj(lP), T(3S) — > 7Xbj(lP), and T(3S) — > IXbj (2P). When these are combined with 
branching ratios for exclusive transitions where the photons from XbJ 7^(15') and 
X6j(2P) — > jT(1S, 2S) and the subsequent decays T(1S, 2S) — > also are observed, one 
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TABLE XV Predicted (jKwong and Rosnerl . 119881 ) and measured (|Yao et all 120061 ) branching ra- 
tios for Xbj(2P) = 2 3 Pj radiative El decays. 



Final Predicted B(%) Measured B (%) 
Level state fKwong and Rosner . 1988 ) fYao et a/. . ,2006) 

2 3 P 7 + 15" 0.96 0.9 ± 0.6 

7 + 25 1.27 4.6 ±2.1 

2 3 Pi 7 + 15 11.8 8.5 ±1.3 

7 ±25 20.2 21 ±4 

2 3 P 2 7 + 15 5.3 7.1 ±1.0 

7 + 25 18.9 16.2 + 2.4 



can obtain branching ratios for the radiative El decays of the Xbj(lP) and -\ h j{2P) states. 



'he Y hj(lP) branching ratios have not changed since the treatment of iKwong and Rosner 



( 119881 ). and are consistent with the predictions quoted there. There has been some improve- 



ment in knowledge of the Xbj{2P) branching ratios, as summarized in Table IXVl 

The dipole matrix elements for T(25) — > 7X W (1P) and T(35) — ► 7x fe j(2P) are shown 
in Figs. [15] and [T6], along with predictions of various models. The dipole matrix element 
predictions are in generally good agreement with the observed values. 

As already pointed out, the most notable exceptions are the matrix elements 



(3 3 5i|r|l 3 Pj). In the NR limit this overlap is less than 5 % of any other S — P over 



suppression occurs for a broad range of potential shapes (IGrant and Rosnerl . Il992l ). This dy 



ap, and its 



namical accident makes these transition rates very sensitive to the details of wave functions 
and relativistic corrections which are not known to this level of precision. This sensitivity is 
shown most clearly looking at the signs of the matrix elements as well as their magnitudes. 
The ave rage experimental v alue for this matrix elemen t is (3 3 Si | r | l 3 Pj) = 0.05 ± 0.006 
GeV -1 ( iCinabro et all 120021 ). Taking the predictions of iGodfrey and Isgurl (119851 ) for com- 
parison, the average over J values gives 0.052 GeV -1 which is in good agreement with the 
observed value. However, more detailed scrutiny gives 0.097, 0.045, and -0.015 GeV -1 for 
J = 2, 1, and matrix elements respectively. Not only is there a large variation in the mag- 
nitudes but the sign also changes, highlighting how sensitive the results for this particular 
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transition are to details of the model due to delicate cancellations in the integral. 

The branching ratios can also be used to measure the ratios of vario us El matrix elements 
which can then be compared to potential model predictions. CLEO (ICinabro et all . 120021 ) 
obtained the following values for ratios: 

0.105 ±0.004 ±0.006, 
0.087 ±0.002 ±0.005, 
0.096 ± 0.002 ± 0.005, 



|(2 3 P 2 


r\l 3 S 1 )\ 


|(2 3 P 2 


r|2 3 £i)| 


|(2 3 P! 


r|l 3 ^i)| 


I<2 3 P 


r|2 3 5i)| 


l<2 3 Pi,2 


r|l 3 ^i)| 


l(2 3 Pi, 2 


r|2 3 5i)| 



where the final ratio averages the results for J = 1 and J = 2. In nonrelativistic calculations 
the El matrix elements do not depend on J. The deviation of the results for J = 1 and 
J = 2 from each other suggests relativistic contributions to the matrix elements. 

D. D-wave states 



The precise information on the masses of 5*-wave and P-wave bb levels leads to 



high 



y constrained predic t ions for the masses and p r oduct ion rates for the P-wave lev- 



els ([Godfrey and Rosner 



2001a 



Kwong and Rosner 



1988). The CLEO Collaboration 



(jBonvicini et all \2004 ) has presented evidence for at least one of these levels in the four- 
photon cascade T(35) -> 7X6(2P), Xb(2P) - lT(lD), T(LD) - W lP), y h (lP) 



yT(l S), followed by the T(IS') annihilation into e + e~ or fi + n~ . CLEO III ( IBonvicini et al. 
20041 ) finds their data are dominated by the production of one T(1P) state consistent with 
the J = 2 assignment and a mass (10161.1 ± 0.6 ± 1.6) MeV, which is consistent with pre- 
dictions from potential models and lattice QCD calculations. The signal product branching 
ratio obtained is £>(7777£ + £ - )y(id) = (2.5 ±0.5 ±0.5) • 10 -5 where the first error is statistical 
and the second one is systematic. The branching ratio is consistent with th e theoretical esti- 



mate of 2.6 x 10 5 ( 1 Godfrey and Rosnerl . 
intermediate state. 



2001a 



Kwong and Rosnerl . Il9881 ) for the T(1 3 P 2 ) 
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FIG. 15 El dipole transition matrix elements for the bottomonium decays 2 3 Si — ► 1 3 -Pj. The 



l abels are the same as in Fig. [3] with t 



(jKwong and Rosner 



le addition of two sets of predictions: KR-Kwong Rosner 



19881 ). F-Fulcher (jFulcherl . 



1990) • 
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E. New hadronic transitions 



1- X&i,2(2P) o/T(15) 

The first transition of one heavy quark onium state to anot her involving uj emission 

2004^ 1: T(2 3 P 1>2 ) -> uT(lS), 



was reported by the CLEO Collaboration (jSeverini et al 



which we have already mentioned in connection with the corresponding transition for the 
Xci(2P) (2 3 Pi) charmonium state. 



2. X61,2(2P) -> X M,2 

The transiti o ns Y fe(2P) — > x&(1P)7T7t have been observed for the first time 
f lCawlfield et all \200& \. One looks for T(35) -> 7X&(2P) -> 77T7tx6(1P) -> 77r7r7T(lS) 



in CLEO data consisting of 5.8 million 3S events. Both charged and neutral pions are de- 
tected. Assuming that r(x w (2P) — > 7T7tx 61 (1P)) = r(% fc2 (2P) — > 7T7tx&2(1P)), both are 
found equal to (0.83 ± 0.22 ± 0.08 ± 0.19) keV, with the uncertainties being statistical, in- 
ternal CLEO systematics, and common systematics from outside sou r ces. T his value is in 



satisfactory agreement with theoretical expectations (IKuang and Yard . 



198J ). 



3. Searches for T(2S, 3S) t/T(1S) 



The decay ip{2S) — > rjJ/ifj(lS) has been known to occur since the early decays of char- 
monium spectro scopy. The world average for its branching ratio is B[ip(2S) — > r)J /ip(lS)\ = 



(3.09 ± 0.08)% (jYao a/.l . 120061 ). The correspon ding T(2S) — » ^Tf lff) process is repre- 



sented by the published upper limit i3 < 2 x 10 3 (iFonseca et al. 
upper limit for T(35) -> r/T(15) is B < 2.2 x 10" 3 (Brock a/ 



984T1 . The corresponding 



19911 ). However, because 



these transitions involve a quark spin flip, they are expected to be highly suppressed in the 
bb system. Defining the ratios 

r[T(2£) - V T(1S)] T[T(3S)^ V T(1S)] 



R' 



r[i>(2S) - »7J/V(15)] 



r[^(25)-^^Mi5)] 



(34) 



Yan (jYanl . Il980l ) estimates P' 
P' = 0.0025, P" = 0.0013. 



1/400, while Kuang (jKuand . 120061 ) finds in one model 
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Combining these results with the latest total widths (jYao et all 120061 ). one predicts 



B[T{2S) -> 7jT(lS)] = (8.1 ±0.8) x 
B[T(3S) -> rjT(lS)] = (6.7 ± 0.7) x 10~ 4 



(35) 
(36) 



The present CLEO III samples of 9 million T(2S) and 6 mi 



used to test these predictions. Preliminary results (IKreinickl . 120071 ) indicate B[T(2S) 



l ion T (3S) decays are being 



v r(lS)\ = (2.5 ± 0.7 ± 0.5) x 10- 4 and B[T(2S) -> 7T°T(1£)] < 2.1 x 10~ 4 (90% c.l.) 



F. Searches for spin-singlets 

Decays of the T(1S', 2S, 3S) states can yield bb spin-singlets, but none has been seen yet. 
One expects IS, 2S, and 3S hyperfine splittings to be approximately 60, 30, 20 MeV/c 2 
([Godfrey and Rosnerl . 12001b). The lowest P-wave singl et state ( u h b ") is expected to be near 
(m(l 3 Pj)) ~ 9900 MeV/c 2 flGodfrev and Rosneil . hooj l. 



The CDF Collaboration has identified events of the form B c — > J/i/jti^, allowing 



Aaltonen et al. 



or the 



fl2007f ). 



first time a precise determination of the mass. The value quoted in 
m(i? n ) = (6275.6±2.9± 2.5) MeV/c 2 , is in reasonable accord with the latest lattice prediction 
of 6304±12+ 48 MeV ( jAUison et all 120051 ). 

The mass of the observed be state can be used to distinguish among various theoretical 
approaches to cc, be, and bb spectra. In this manner, in principle, one can obtain a more 
reliable prediction of the masses of unseen bb states such as r)b(lS, 2S, 3S). For example, 
by comparing predictions of potential models to the measured values of the J/t/> , r? r , T 



and B r states one could use the prediction o: 



Eichten and Quigg 



1994 



Fulcher 



1991 



1999 



the most reliable models ( lEbert et al. 



Godfrey and Isgurl . ll985l ) to estimate the mass 



2003a ; 



of the r) b (lS) = 9400 - 9410 MeV. 

Several searches have been performed or are under way in IS, 2S, and 3S CLEO data. 
The allowed Ml transition in T(IS') — > r yr]b(lS) can be studied by reconstructing exclusive 
final states in rjb(lS) decays. One may be able to dispense with the soft photon, which could 
be swallowed up in background. Final states are likely to be of high multiplicity. 

One can search for higher-energy but suppressed Ml photons in T(n'S) — > / -fi]i ) (nS) 
(n ^ n') decays. Inclusive searches already exclude many m odels. The strongest upper limit 
obtained is for n' = 3, n = 1: B < 4.3 x 10~ 4 (90% c.l.) (jArtuso et all |2005|). Exclusive 
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searches (in which rjb decay products are reconstructed) also hold some promise. Searches 
for rjb using the sequential processes T(3S) — + ^l hh (l 1 P^ ) — > 7i°'jrjb(lS) and T(35) — > 
7Xbo(2-P) — > r yrjrjb(lS) (suggested in IVoloshinl (12004 )) are being conducted. Additional 
searches for hb involve the transition T( 3S) — > 7r + 7r~ hb [for which a typical e xperimental 



1991 



Butler et al. 



upper bound based on earlier CLEO data (IBrock et al. 
The hb — > 7?7b transition is expected to have a 40% branching ratio ( 
20021 ) . much like h c — ► 7?7 C . 



1994h is Q(1Q- 3 )] 



Godfrey and Rosnerl . 



G. T(4S) 



The T(4S) is the lowest-lying bound bottomonium state above open-flavor threshold. 



Its mass and total width as wel 



recently by BaBar (jAubert et al. 



as elec tronic width have been determined in scans, most 



2005ah : M = (10579.3±0.4±1.2) MeV/c 2 , T ee = (0.321 ± 



0.017 ± 0.029) keV, T = (20.7 ± 1.6 ± 2.5) MeV. Although the T(45) has primarily been 
regarded as a BB "factory," its decays to bound hb states are beginning to be observed 
in the large data samples accumulated by BaBar and Belle. This is not surprising, as the 
corresponding first char monium state ab ove flavor threshold, the ^(3770), does decay - 



rarely - to charmonium (lYao et al. 



2006). 



The BaBar Collaboration (jAubert et all l2006bl ) measures the product branching frac- 
tions 

B[T(4S) -> 7i + 7i-T(lS)] x B(T(15) -> fj+fi-) = (2.23 ± 0.25 ± 0.27) x lO" 6 and 



B[T(AS) -> 7T + 7r-T( 26')1 x B(T(2S ) 



ji ji 

fJL + (JL-) -- 



1.69 ± 0.26 ± 0.20) x 10~ 6 , while the 



20071) finds B[T(4S) -> n+n-T(lS)] x B(T(15') 



Belle Collaboration (jSokolov et al. 
yU + /i~) = (4.4 ± 0.8 ± 0.6) x 10^ 6 . These product branching fra ctions, when co mbined 
with B{T(1S)[T{2S)} -> ji + ji-) = (2.48 ± 0.05)%[(1.93 ± 0.17)%] JYao et all \200& \ result 
in branching fractions of the order of 10~ 4 and partial widths of a few keV, comparable 
with other partial widths for dipion transitions in the T system of the same order of 
magnitude. An interesting feature is that the distribution of m(7r + 7r~) in T(4S f ) — > T(2S) 
looks markedly different from the T dipion transitions with An = 1 [T(35) — > T(2S*), 
T(2S) -> T(1S)] and more resembles that of T(35) -> T(1S); however, the T(4S) -> T(1S) 
dipion spectrum (An = 3) can be described by a mod el that suits the A n = 1 bottomonium 



transitions and also the shape in ip(2S) — ► 7r + 7r J/tp (jKuang and Yanl . Il98ll ) 
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FIG. 17 Invariant mass of th e dipion sys tem in Y(4S) 



BaBar ( Aubert et al. 



20061)1 ) and Belle (|Sokolov et al 



7T + 7 r T(1,2S') as measured in data from 



20071 ) (points), after efficiency correction. 



F or the BaBar figu r es, th e dotted line is the selection efficiency, and the solid line is the prediction 



of 



Kuang and Yanl (|198ll ) . In the Belle plo t, the shaded histogram i s a background e stimate, an d 



the curve is based on the model detailed in 



Brown and Cahnl (|1975l ); 



Voloshin 



(|1975T ): 



Yanl dl980). 



The measured dipion invariant mass distributions for T(45) — > 7t + tt T(1S, 25) are shown 
in Figure [T71 



H. States above open flavor threshold 



Two states have been seen in e + e scattering (jYao et all 120061 ). establishing quantum 
numbers J PC = 1 — ; T(10860) (mass 10.865 ± 0.008 GeV, total width 110 ± 13 MeV) 
and T(11020) (mass 11.019 ± 0.008 GeV, total width 79 ± 16 MeV). These states are often 
identified as 5S and 6S bottomonium levels. 



VI. SUMMARY 

In the presence of much more accurate data, multipole expansions for both electromag- 
netic and hadronic transitions hold up well. The coefficients appearing in these expansions 
have been described in the past by a combination of potential models and perturbative 
QCD. As expected there are significant relativistic corrections for the charmonium system. 
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The overall scales of these corrections are reduced for the bb system and are consistent with 
expectations from the NRQCD velocity expansion. Relativistic corrections are determined 
in the same framework as leading order terms. However, relativistic corrections have not 
improved markedly upon the nonrelativistic treatments, though some qualitiative patterns 
(such as hierarchies in electric dipole matrix elements) are reproduced. 

Electromagnetic transitions for which the leading-order expansion coefficient is dynam- 
ically suppressed are particularly sensitive to relativistic corrections. For the T(3S) — > 
XbO-P) El transitions there is a large cancellation in overlap amplitude because of the node 
in the 3S* radial wavefunction. The result is a wide scatter of theoretical predictions. For 
the T(3S) — > rj^lS) Ml transition, the overlap coefficient vanishes in leading order (a hin- 
dered transition). Here the experimental upper bound on the rate is smaller than expected 
in potential models for relativistic corrections. Modern theoretical tools (effective theories 
and nonperturbative lattice QCD) combined with more detailed high-statistics experimental 
data will help pin down the various relativistic corrections. 

Decays described by perturbative QED or QCD, such as Xco,2 — > (77, gg), appear to 
behave as expected, yielding values of as for the most part consistent with other determi- 
nations. Exceptions (as in the case of the anomalously small J/ip hadronic width) can be 
ascribed to large QCD or relativistic corrections or to neglected color-octet components of 
the wave function which are not yet fully under control. 

Recent experiments have also observed a number of new hadronic transitions. Many 
details remain to be understood. The two-pion invariant mass distributions in both the 
T(3S) — > T(IS') + 2n and T(AS) — > T(2S') + 2n transitions do not show typically strong 
S'-wave behavior. Perhaps some dynamical suppression plays a role in these transitions. To 
further complicate the situation, the T(4S) — > T(1S) + 2tt decay seems to show the usual 
S'-wave behavior with the dipion spectrum peaked toward the highest effective masses. 

Coupled-channel effects appear to be important in understanding quarkonium behavior, 
especially in such cases as the X(3872) which lies right at the D°D*° threshold. It seems that 
long-awaited states such as "molecular charmonium" [with X(3872) the leading candidate] 
and hybrids [perhaps such as Y(4260)] are making their appearance, and the study of their 
transitions will shed much light on their nature. Now that we are entering the era of precise 
lattice QCD predictions for low-lying quarkonium states, it is time for lattice theorists to 
grapple with these issues as well. 
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